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Recent  advances  in  sensor  and  guidance  system  technology  coupled  with 
projected  improvements  in  cluster  munitions  apoear  to  offer  the  possibility  of 
delivering  effective  nonnuclear  payloads  to  within  a few  meters  of  targets  that 
are  located  hundreds  of  kilometers  from  the  initial  launch  point.  Using  unmanned 
delivery  systems  in  high  attrition  environments  is  particularly  attractive  since 
effective  munitions  delivery  against  high  value,  preprogrammed  targets  can  be 
accomplished  without  exposing  expensive  aircraft  to  severe  enemy  air  defenses. 

Because  of  the  attractiveness  of  unmanned,  deep  strike  missions,  many  standoff 
missile  (SOM)  and  remotely  piloted  vehicle  (RPV)  concepts  have  been  proposed  by 
various  OOD  organizations  and  defense  contractors.  The  Strike  Options  Comparison 
(SOC)  Study  was  initiated  by  HQ  AFSC/XR  to  provide  insights  into  the  develooment 
planning  needs  associated  with  long  range  munitions  delivery.  The  SOC  Study  was 
structured  with  the  intent  of  providing  consistent  and  realistic  cost/effective 
comparisons  of  various  SOM  and  RPV  concepts  for  accurately  delivering  munitions 
on  targets  as  far  as  650  kilometers  beyond  the  forward-edge-of-the-battl e-area 
(FEBA).  As  a benchmark,  the  capabilities  of  the  F-lllD  in  performing  the  deep 
strike  mission  were  subjected  to  the  same  comparison  analysis  as  the  SOM  and  RPV 
conceptual  systems. 

The  diverse  characteristics  of  the  strike  options  considered  in  the  study 
presented  a broad  spectrum  of  survivability  problems.  One  of  the  key  analytical 
tasks  in  the  SOC  Study  was  to  obtain  consistent  survivability  estimates  for  the 
various  deep  strike  delivery  systems  penetrating  areas  defended  by  surface-to-air 
missiles  (SAM),  airborne  interceptors,  and  terminal  defenses.  In  this  report,  a 
simplified  methodology  for  estimating  survivability  of  deep  strike  raids  penetrat- 
ing through  SAM  area  defenses  is  described.  A listing  of  a digital  computer  model 
which  was  developed  to  implement  the  methodology  is  provided  in  an  appendix  to 
the  report.  The  material  presented  herein  was  carefully  selected  to  avoid  security 
classification  and  thus  allow  the  methodology  and  model  oescriotion  to  be  dis- 
tributed without  undue  hindrance.  Classified  data  relating  to  survivability  of 
SAM  area  defenses,  as  well  as  airborne  interceptors,  and  terminal  defenses  are 
contained  in  DAS-TR-77-8,  Strike  Options  Comparison  (SOC)  Study  - Final  Report, 
SECRET,  publication  pending. 
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SECTION  1 


INTRODUCTION 

The  probability  of  an  airborne  vehicle  surviving  a penetration  of  territory 
defended  by  hostile  ground  based  weapons  is  dependent  upon  many  factors.  These 
factors  may  generally  be  thought  of  as  influencing  either  the  number  of  defenses 
encountered  or  the  outcome  of  an  encounter  when  it  takes  place.  The  analyst  has 
the  choice  of  modeling  these  factors  either  deterministically  or  stochastically. 

In  the  stochastic  approach,  the  problem  is  usually  modeled  in  elaborate  detail, 
with  the  end  result  being  the  average  of  many  trials  (replications)  based  on 
"MONTE  CARLO”  random  draws.  A stochastic  model  is  generally  justified  only  if 
the  nature  of  the  statistical  distributions  describing  the  modeled  random  processes 
can  be  estimated  and  are  believed  to  be  important  to  the  final  answer.  A determin- 
istic (expected  value)  model,  on  the  other  hand,  is  based  only  on  the  statistical 
means  of  the  distributions  describing  the  random  model  processes.  As  a rule,  expec- 
ted value  models  require  fewer  and  less  detailed  inputs  than  a stochastic  model, 
and  they  produce  many  more  results  than  a stochastic  model  per  unit  analyst  time 
and  per  unit  computer  time.  For  these  reasons,  the  deterministic  approach  is  by 
far  the  superior  of  the  two  for  situations  involving  inadequately  defined  systems 
or  processes.  Such  situations  are  overwhelmingly  pervasive  in  the  authors' 
experience. 

The  FIRE  survivability  model,  developed  for  use  in  the  Strike  Options  Compari- 
son (SOC)  Study  by  the  Directorate  of  Aerospace  Studies,  is  a deterministic  model. 
Its  approach  is  straightforward:  (1)  determine,  for  the  given  conditions,  the 
average  probability  of  a penetrator  in  a raid  surviving  an  encounter  with  a single 
site  of  each  type  of  SAM  (the  average  being  taken  with  respect  to  penetrator  off- 
set relative  to  the  site),  then  (2)  determine  the  subsequent  expected  survivability 
of  a member  of  the  raid  as  a function  of  penetration  distance  (sites  potentially 
encountered). 

The  model  contains  provisions  for  handling  in  varying  detail  the  following 
aspects  of  the  survivability  problem  relative  to  ground  based  defenses: 


3 


Penetration  Scenario 

Corridor  dimensions  (the  area  containing  the  defenses) 
Penetration  distance 

Distribution  of  defenses  (up  to  10  types  of  SAMs) 

Raids  of  Penetrators 

Number  of  penetrators 

Two-dimensional  distribution  of  penetrators 
Penetration  altitude 
Penetration  speed 

Penetration  Aid  Effectiveness  (heuristic) 

Decoys 

Chaff 

Electronic  Countermeasures  (ECM) 


Defense  Availability 

Terrain  degradation  of  firing  opportunities 

Reduction  of  number  of  defenders  from  previous  suppression  efforts 
Defense  readiness  (to  exclude  travel  time,  downtime,  etc.) 

Surface-to-Air  Missile  (SAM)  Characteristics 
Launchers  per  site 
Missiles  per  launcher 
Launcher  reload  time 
Radar  self-screening  angle 
Maximum  effective  flyout  range 
Average  missile  flyout  speed 
Probability  of  kill 
Sensor  height  above  ground 
Maximum  sensor  range 
Geometric  missile  launch  restrictions 
Time  from  raid  detection  to  first  salvo 
Time  interval  between  salvos 
Salvo  reliability 


The  manner  in  which  the  previously  mentioned  factors  are  considered  in  the 
model  will  be  discussed  in  section  2 with  references  to  the  associated  input 
variables.  The  input  technique  and  the  input  variables  are  discussed  in  section 

3.  The  calculation  of  the  penetrator  survival  probability  is  treated  in  section 

4,  while  the  calculation  of  the  number  of  salvos  fired  at  the  raid  for  a given 
offset  is  explained  in  section  5.  Two  appendices  conclude  the  report.  Appendix 
A discusses  the  difficulties  of  modeling  the  effects  of  terrain  masking  relative 
to  low  altitude  penetrators.  Appendix  B provides  a well -commented  FORTRAN  listing 
of  the  model. 
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SECTION  2 
MODEL  DESCRIPTION 


A.  PENETRATION  SCENARIO 

The  FIRE  model  uses  the  concept  of  a rectangular  penetration  corridor 
(CORDWTH,  NDIST)^  in  which  all  SAM/penetrator  interactions  take  place.  Nominally 
the  model  corridor  is  associated  with  some  real  geographic  area.  The  number, 
distribution,  and  capabilities  of  SAM  defenses  within  the  model  corridor  are 
based  on  an  estimate  of  what  is  expected  to  be  in  the  geographic  area  in  the  time 
frame  of  the  study.  In  conjunction  with  the  corridor  concept  a partitioned, 
discrete  distribution  of  defenses  is  assumed.  The  area  of  the  corridor  behind 
the  FEBA  is  divided  into  adjacent  intervals  each  n kilometers  deep  (DEL),  and  the 
number  of  each  type  of  SAM  is  specified  for  each  iriterval  (NSAML,  NSAMH).  Within 
a given  interval  locations  are  not  specified  but  the  SAMs  are  assumed  to  be 
randomly  positioned  in  a direction  parallel  to  the  FEBA. 

Two  distinct  sets  of  up  to  10  SAM  types  may  be  defined  in  the  model  at  any 
one  time.  Usually  one  set  is  defined  as  the  "lev;  altitude"  defense;  the  other 
as  the  "high  altitude"  defense  in  the  corridor.  This  aspect  of  the  program  is 
merely  a convenience,  for  each  set  is  identically  treated  in  the  program  logic. 
Both  sets  could  as  well  represent  low  or  high  altitude  defenses  (I LOW). 

B.  PENETRATOR  RAIDS 

Penetrators  may  fly  singly  or  in  groups  (raids)  of  two  or  nore  (NPENT).  The 
FIRE  model  allows  only  a single  penetrator  type  to  be  considered  in  a given  case. 
The  survivability  of  a penetrator  is  very  much  dependent  upon  the  size  of  the 


1.  Throughout  this  section,  references  to  input  identifiers  are  given  to  direct 
the  reader  to  the  appropriate  input  cards  defining  the  aspect  of  the  model  under 
discussion.  These  identifiers  and  associated  inputs  are  discussed  in  the  follow- 
ing section.  Not  all  Identifiers  are  discussed  in  this  section.  The  reader's 
attention  is  specifically  directed  to  the  "bookkeeping"  Identifiers  TITLE,  IPINT, 
ENDCASE,  and  ENOJOB  as  well  as  NPHI,  TERINO,  and  XLAMBDA  which  are  discussed  in 
section  4 or  5. 
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group  of  which  it  is  a part  because  of  the  dilution  of  the  defenses.  For  example, 
a single  penetrator  receives  all  the  fire  from  a SAM  site;  a penetrator  in  a 
group  of  10  receives  on  the  average  about  1/10  of  the  fire.  A raid  may  be  started 
on  the  enemy  side  of  the  FEBA  to  simulate  a raid  of  standoff  weapons  launched 
behind  the  FEBA  (PSTART). 

A raid  composed  of  multiple  penetrators  must  have  dimension.  Ideally,  it 
should  be  structured  in  three  dimensions.  The  FIRE  model  allows  only  for  raid 
width  (SP)  perpendicular  to  the  direction  of  flight  and  raid  depth  (TR)  parallel 
to  the  direction  of  flight.  The  raid  is  assumed  to  be  rectangular  and  to  be  con- 
tained in  a plane  parallel  to  the  ground.  The  raid  width  affects  the  probability 
of  encountering  a SAM  site;  the  depth  affects  the  length  of  time  the  raid  is 
exposed  to  fire  from  a site.  Within  the  model,  the  dimensions  of  the  rectangle 
are  determined  by  the  number  of  penetrators  and  the  number  of  equally  spaced 
columns  and  rows  into  which  they  are  placed  (PENC).  The  situation  is  depicted  in 
figure  1.  The  penetrator  speed  (VP)  affects  the  length  of  time  a raid  is  exposed 


NUMBER  OF 
COLUMNS  AND 
SPACING  (INPUT) 

NUMBER  OF  ROWS  (CALCULATED) 
SPACING  OF  ROWS  (INPUT) 


Figure  1.  Raid  Structure 


to  fire.  This  In  turn  affects  the  number  of  salvos  which  may  be  fired;  hence,  the 
probability  of  penetrator  survival.  Penetrator  altitude  (HM)  also  affects  the  time 
a raid  is  exposed  to  fire  from  a site.  Reducing  penetrator  altitude  has  the  effect 
of  postponing  the  point  of  detection  (see  section  5);  hence,  reducing  the  exposure 
time  of  the  penetrator  to  the  defenses.  The  effect  is  more  pronounced  the  lower 
the  penetrator. 
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C.  PENETRATION  AIDS 


Penetration  aids  can  reduce  the  number  of  salvos  fired  and/or  the  effectiveness 
of  salvos  which  are  fired  (PKH,  PKL). 

The  accurate  determination  of  the  effectiveness  of  penetration  aids  is  nearly 
an  impossible  task.  Furthermore,  in  most  comparative  or  parametric  evaluations  of 
survivability,  the  fundamental  question  concerns  what  may  be  gained  from  a penetra- 
tion aid  of  some  preselected  effectiveness.  The  approach  taken  in  FIRE  is  simply 
to  allow  the  specification  of  effectiveness  of  the  penetration  aid  over  some 
selectable  portion  of  the  penetration.  For  example,  it  is  possible  to  specify  that 
chaff  is  used  (ICHAFF)  over  some  interval  of  the  penetration  (CHAFFS,  CHAFFE),  and 
that  it  is  252  effective  in  that  interval  (CHEFF)  in  that  it  reduces  the  number  of 
effective  SAM  encounters  to  75%  of  what  they  would  be  without  chaff.  Chaff  effective- 
ness can  be  specified  by  SAM  type,  but  the  one  interval  which  may  be  specified  applies 
to  all  SAM  types. 

In  the  case  of  ECM  (lECM),  the  effects  are  treated  through  a reduction  of  the 

SAM  probability  of  kill  (ECMEFF).^  For  example,  if  ECM  is  25%  effective,  it  reduces 

the  single  shot  kill  probability  to  75%  of  its  original  value.  As  with  chaff,  ECM 

may  be  specified  over  some  fraction  of  the  penetration  distance  (ECMS,  ECME).  ECM 

effectiveness  also  can  be  specified  by  SAM  type,  and  the  one  interval  which  may  be 
specified  applies  to  all  SAM  types. 

The  final  penetration  aid  which  is  modeled  is  the  decoy  (IDECOY).  Decoys  are 
assumed  to  be  perfect  in  the  sense  that  the  defense  is  as  likely  to  fire  at  a 
decoy  as  at  a lethal  penetrator.  Thus,  their  effect  is  to  dilute  the  defenses  as 
would  additional  penetrators.  Decoys  may  be  specified  over  some  fraction  of  the 
penetration  (DSTR,  DENO).  They  do  not  enter  into  the  calculation  of  the  raid 
dimensions  in  the  model.  If  it  is  desired  to  treat  the  decoys  as  expanding  the 


2.  Deceptive  ECM  will  lower  the  P,  of  a missile.  Noise  ECM  will  behave  in  a manner  * 
similar  to  chaff  by  reducing  the  nOmber  of  effective  SAM  encounters  or  shortening  ^ 
the  encounter  by  delaying  burnthrough,  thus  reducing  the  number  of  salvos  fired.  | 
This  latter  effect  may  be  simulated  by  reducing  the  radar  range  or  effective  SAM  | 
range  as  deemed  appropriate.  ; 
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raid  dimensions,  the  spacing  given  for  the  lethal  penetrators  may  be  increased 
aporopriately.  Only  one  type  of  decoy  may  be  considered,  although  it  may  be  used 
to  decoy  multiple  SAM  types  as  specified  (PENCDEC,  OECPPEN).  Penetrators 
carrying  decoys  are  attrited  in  proportion  to  their  number. 

D.  SAM  CHARACTERISTICS 

The  number  of  missile  salvos  which  a SAM  site  can  fire  at  a raid  is  deter- 
mined primarily  by  the  time  the  penetrators  are  exposed,  the  number  of  salvos 
available  for  firing,  and  the  rapidity  with  which  they  can  be  launched.  FIRE 
permits  specifying  many  parameters  affecting  these  areas. 

The  number  of  missiles  per  launcher  (NML)  and  the  number  of  launchers  per 
site  (ISET,  ISETL)  are  specified  as  well  as  the  number  of  missiles  constituting 
a salvo  (SALSZ).  Additionally,  the  time  from  penetrator  detection  to  firing  the 
initial  salvo,  the  subsequent  time  between  salvos  (ISET,  ISETL),  and  the  launcher 
reload  time  (SRLD)  are  specified.  Clearly  the  flyout  time  of  the  missile  salvo 
is  also  important.  Missile  flyout  time  is  calculated  in  FIRE  on  the  basis  that 
the  missile  always  flies  to  the  intercept  point  in  a straight  line  at  constant 
speed  (ISET,  ISETL).  Although  this  is  an  unsophisticated  model,  it  is  generally 
adequate.  Suppose,  for  example,  a high  altitude  intercept  takes  place  at  20  km 
from  the  site  and  that  the  actual  flyout  time  is  40  sec  (true  average  straight 
line  speed  of  500  m/sec).  If  the  average  estimated  speed  used  in  the  model  were 
600  m/sec,  the  model  flyout  time  would  be  about  33  seconds  — a 7 second  error. 
Such  an  error  is  relatively  unimportant  when  considered  with  uncertainties  in 
firing  doctrine,  time  between  salvos,  terrain  masking,  site  availability,  etc. 

The  SAM  sensor  as  modeled  in  FIRE  is  specified  by  three  parameters  --  sensor 
height  above  ground,  maximum  sensor  range  (ISET,  ISETL),  and  sensor  self-screening 
angle  (ALPHA).  The  self-screening  angle  describes  how  close  to  the  terrain  the 


3.  The  time  between  salvos  is  the  time  from  "Intercept"  by  one  salvo  to  the 
launching  of  the  next.  When  multiple  missile  salvos  are  fired,  there  is  an 
Interval  between  the  individual  missile  launches.  FIRE  does  not  explicitly  model 
this  interval,  but  it  may  be  allowed  for  by  adjusting  the  time  between  salvos. 
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radar  may  look.  It  may  also  be  used  to  impose  a constant  masking  angle  about  the 
sensor.  For  low  altitude  penetrators,  this  angle  usually  determines  the  radar 
detection  range;  hence,  the  exposure  duration  of  a penetrator  to  the  SAM.  The 
geometry  of  the  situation  is  discussed  and  illustrated  in  section  5.  The  effect 
of  sensor  iieight  is  alsr  treated  in  section  5. 


Initially  terrain  was  brought  into  the  problem  through  a probability  that 
terrain  masking  does  not  prevent  a SAM  from  firing  (what  might  be  called  the  "prob- 
ability of  favorable  terrain")  (ISET,  ISETL).  The  probabilities  were  taken  from 

4 

work  done  for  the  SOC  Study  by  Dr.  Charles  C.  Carson  of  Sandia  Laboratories.  Later 
consideration  of  the  terrain  problem  resulted  in  an  improved  method  for  treating 
terrain  with  the  FIRE  model. 

The  new  methodology  uses  a weighted  average  masking  angle  and  sets  to  unity 

c 

P^,  the  probability  of  favorable  terrain.  This  implies  the  availability  of  a 
probability  distribution  of  masking  angles,  but  calculation  of  meaningful  P^'s 
has  the  same  implication.  The  calculation  of  a weighted  average  masking  angle 
may  be  made  in  several  ways.  The  most  straightforward  method  simply  uses  the 
probabilities  of  occurrence  of  various  masking  angles  as  the  weights.  Another  way 
involves  the  computation  of  the  expected  number  of  shots  per  site  as  a function  of 
constant  masking  angle  and  the  subsequent  determination  of  the  masking  angle  corres- 
ponding to  the  weighted  average  of  expected  number  of  shots. 

Each  SAM  type  in  FIRE  is  modeled  as  being  able  to  fire  at  penetrators  within 
a cylinder  (ISET,  ISETL)  whose  axis  passes  through  the  site  perpendicular  to  the 
earth.  A smaller  coaxial  cylinder  defines  a "dead  zone"  (ISET,  ISETL)  within  which 
penetrators  are  safe  from  intercept.  The  dead  zone  arises  primarily  from  launcher 
elevation  limitations  and  warhead  fuzing  delays. 


f-i 


4.  Taken  from  a classified  report. 

5.  During  the  course  of  the  SOC  Study,  the  only  masking  angle  used  was  the  sensor 
self-masking  angle. 
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The  model  allows  the  first  salvo  to  be  fired  a specified  time  after  the 
penetrator  enters  site  sensor  coverage  (ISET,  ISETL).  Sensor  coverage  for  radar 
guided  SAMs  is  modeled  as  a hemisphere.  It  can  be  modeled  differently  for  IR  SAMs 
to  account  for  the  fact  that  much  of  the  IR  signal  may  be  masked  as  a penetrator 
approaches  the  site.  This  has  been  handled  by  allowing  detection  to  occur  only 
when  the  azimuth  angle  of  the  site  as  seen  from  the  penetration  reaches  some 
predetermined  value  between  0 and  90  degrees  (ASPECT).  The  situation  is  depicted 
in  figure  2. 

Radar  guided  SAMs  are  restricted  in  the  model  from  launching  at  the  penetrator 
after  the  penetrator  has  passed  the  site.  The  IR  SAMs  are  not  subject  to  this 
restriction.  (The  differentiation  between  radar  guidance  and  IR  guidance  was  made 
after  the  SOC  Study.) 


Figure  2.  The  Definition  of  the  Input  Parameter  ASPECT 
E.  AVAILABILITY  OF  DEFENSES 

Defenses  in  the  field  are  available  only  a percentage  of  the  time.  They  are 
subject  to  downtime  for  mechanical  and  electronic  problems,  they  may  be  in  transit, 
or  their  crews  may  not  be  capable  of  responding  in  time  for  various  reasons.  The 
FIRE  model  allows  for  a "readiness  factor"  which  is  the  probability  a defense 
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site  is  available  (ISET,  ISETL).®  It  is  also  possible  to  specify  the  probability 
that  a launcher  at  a site  is  ready  (PLCHRDY)  and  the  reliability  of  a salvo  (ROS). 
Salvo  reliability  is  a multiplier  of  missile  salvo  Pj^. 

In  addition  to  the  reaciness  factors,  FIRE  allows  the  availability  of  each 
SAM  type  to  be  decreased  by  an  arbitrary  factor  (DF).  This  allows  the  analyst, 
for  example,  to  simulate  some  level  of  defense  suppression  for  the  various  SAMs 
without  having  to  explicitly  respecify  their  distribution. 

It  is  important  to  note  that  the  model  has  no  altitude  or  speed  restrictions 
by  which  SAM  types  are  judged  ineligible  to  fire  at  a raid.  SAM  types  which 
cannot  fire  at  a raid  for  either  of  these  reasons  should  not  be  input. 


6.  The  readiness  factors  used  in  the  SOC  Study  were  adopted  from  a classified 
reference. 
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SECTION  3 


INPUT 

The  Input  routine  for  FIRE  is  based  on  a concept  first  adopted  by  the  authors 
in  1965.  This  concept,  which  identifies  individual  input  variables  and  groups  of 
variables  by  an  alphanumeric  identifier,  provides  the  user  with  (1)  a printed 
record  of  all  inputs,  (2)  no  requirement  for  ordering  inputs  (except  when  two  or 
more  cards  are  associated  with  an  identifier),  (3)  variation  of  the  initial  data 
set  by  inputing  only  the  changes  (simplifies  stacking  cases),  and  (4)  a check  of 
all  identifiers  against  a dictionary  of  allowed  identifiers,  thus  helping  to 
minimize  input  errors. 

An  identifier  card  always  contains  the  identifier  left  justified  in  column  1. 
Any  remaining  fields  on  the  card  may  contain  data.  The  field  structure  can  vary 

O 

with  the  program;  FIRE  using  the  format  A7,  3X,  E20.8.  The  identifier  is  checked 
for  validity,  and  if  valid,  it  is  used  as  a signal  to  the  program  describing  what 
information  is  on  the  card  and/or  what  associated  cards  should  follow.  All  cards 
read  are  displayed  in  the  output.  Invalid  identifiers,  should  any  be  found,  are 
indicated  and  the  program  is  halted  aftpr  checking  the  nest  of  the  input  for  that 
case. 


There  follows  a list  of  identifiers  used  in  FIRE.  Their  functions  are  des- 
cribed through  a discussion  of  the  input  parameters  associated  with  them.  Param- 
eter default  values  are  given  when  appropriate. 

ALPHA  Indicates  that  the  i-th  field  of  the  following  card  contains 
the  radar  self-screening  angle  for  SAM  type  i which  specifies 
how  close  to  the  terrain  the  radar  can  operate  (default  0.5 
degrees).  This  angle  can  also  be  used  to  simulate  a constant 
masking  angle.  Following  card  has  a 10F8.0  format. 


8.  A useful  format  in  many  instances  is  AlO,  215,  6F10.0.  The  integer  fields  may 
be  used  for  Integer  input  or  for  indices  of  arrays. 
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ASPECT  Defines  the  forward  hemisphere  angular  tracking  capability  of 
SAM  sensors  (see  figure  2,  section  2).  The  default  is  zero 
which  corresponds  to  the  standard  situation  for  radar. 

CHAFFE  Indicates  the  penetrator  distance  in  kilometers  from  the  FEBA 
at  which  the  chaff  effectiveness  modeling  ends. 

CHAFFS  Indicates  the  penetration  distance  in  kilometers  from  the  FEBA 
at  which  the  chaff  effectiveness  modeling  begins. 

CHEFF  Indicates  that  the  i-th  field  of  the  following  card  contains  the 
chaff  effectiveness  for  SAM  type  i.  A value  of  zero  indicates 
that  the  chaff  is  completely  effective.  A value  of  one  indicates 
no  effectiveness.  This  parameter  is  a multiplier  of  the  proba- 
bility of  encountering  a site.  Format  of  following  card  is  10F8.0. 

CORWDTH  Indicates  that  the  first  field  of  the  next  card  specifies  the 

corridor  width  for  the  "low  altitude"  SAMs,  the  second  field  for 
the  "high  altitude"  SAMs  (format  2F8.0). 

OECPPEN  Indicates  that  the  i-th  field  of  the  folVowing  card  contains  the 
number  of  perfect  decoys,  N,  relative  to  SAM  type  i carried  per 
decoy-carrying  penetrator.  Entries  should  be  either  zero  or  N 
as  only  one  type  decoy  is  permitted.  Following  card  has  a 10F8.0 
format. 

DEL  The  depth  in  kilometers  of  the  intervals  specifying  the  positions 
of  the  defenses  (default  10). 

OENO  Indicates  the  penetration  distance  In  kilometers  from  the  FEBA 
at  which  the  decoy  effectiveness  modeling  ends. 

DF  Indicates  that  the  1-th  field  of  the  next  card  contains  a multi- 

plier of  the  availability  of  SAM  type  1.  The  multiplier  repre- 
sents the  probability  a SAM  site  has  been  suppressed.  The 
format  of  the  following  card  Is  10F8.0.  (Default  1.0.) 
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DSTR  Indicates  the  penetration  distance  in  kilometers  from  the  FEBA 

at  which  the  decoy  effectiveness  modeling  begins. 

ECME  Indicates  the  penetration  distance  in  kilometers  from  the  FEBA 

at  which  the  ECM  effectiveness  modeling  ends. 

ECMEFF  Indicates  that  the  i-th  field  of  the  following  card  contains 

the  ECM  effectiveness  for  SAM  type  i.  A value  of  zero  indicates 
that  the  ECM  is  completely  effective.  A value  of  one  indicates 
no  effectiveness.  This  parameter  is  a multiplier  of  the  single 
shot  P|^  (not  the  salvo  which  is  input).  Format  of  following 
card  is  10F8.0. 

ECMS  Indicates  the  penetration  distance  in  kilometers  from  the  FEBA 

at  which  the  ECM  effectiveness  modeling  begins. 

ENDCASE  Signifies  the  end  of  input  for  the  case  in  question. 

ENDJOB  Signifies  the  end  of  a run.  The  ENDJOB  card  does  not  perform 

the  function  of  ENuCASE  card. 

HM  The  altitude  of  the  penetrator  above  the  ground  measured  in 

kilometers. 

ICHAFF  A value  of  unity  indicates  that  chaff  effectiveness  will  be 

modeled  for  the  present  case.  A value  of  zero  indicates  other- 
wise (default  zero). 

IDECOY  A value  of  unity  indicates  that  decoy  effectiveness  will  be 

modeled  for  the  present  case.  A value  of  zero  indicates  other- 
wise (default  zero). 

lECM  A value  of  unity  Indicates  that  ECM  effectiveness  will  be  modeled 
for  the  present  case.  A value  of  zero  indicates  otherwise 
(default  zero). 
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ILOW  For  a value  of  unity,  SAM  parameters  representing  SAM  charac- 

teristics against  "low  altitude"  penetrators  are  chosen.  For 
a value  of  zero  SAM  characteristics  against  "high  altitude" 
penetrators  are  selected  (default  zero). 

IPINT  A variable  controlling  output.  The  survivability  output  is 
given  for  every  IPINT  times  DEL  kilometers  where  DEL  is  the 
depth  of  the  corridor  intervals  (default  five). 

ISET  Indicates  that  the  following  card  will  specify  for  which  high 

altitude  SAMs  data  are  to  be  read,  and  that  the  cards  following 
that  will  contain  the  SAM  data.  The  format  for  the  card  follow- 
ing the  ISET  card  is  1011.  A "0"  in  column  i indicates  no  data 
for  SAM  type  i are  to  be  read,  a "1"  indicates  that  data  are  to 
be  read.  The  cards  following  this  card  (one  card  for  each  "1" 
and  ordered  by  SAM  type  i)  specify  by  field 

( 1 ) SAM  type  i 

(2)  Number  of  launchers  per  site 

(3)  Probability  that  type  i SAM  site  is  ready  to  fire  if 
called  upon 

(4)  Probability  of  "favorable  terrain" 

(5)  Time  from  first  possible  detection  of  raid  to  firing 
first  salvo  (seconds) 

(6)  Time  between  salvos  (seconds) 

(7)  Maximum  useful  missile  flyout  range  (km) 

(8)  Average  missile  flyout  velocity  (km/sec) 

(9)  Radius  of  dead  zone  (km) 

(10)  Maximum  sensor  range  (km) 

(11)  Height  of  sensor  above  ground  (km) 

The  format  is  12,  13,  3X,  8F8.2,  F8.3. 

ISETL  This  Identifier  and  its  attendant  cards  are  similar  to  ISET  and 
its  associated  cards.  They  describe  SAM  characteristics  against 
low  altitude  penetrators. 


16 


NDIST 

NML 

NPENT 

NPHI 

NSAMH 

NSAML 

PENC 

PENCOEC 


The  number  of  points  spaced  DEL  kilometers  apart  (where  DEL 
is  the  depth  of  the  corridor  intervals)  at  which  survival 
probability  is  to  be  calculated.  NDIST  may  be  regarded  as 
specifying  corridor  depth  or  alternately  maximum  penetration 
distance  (default  65), 

Indicates  that  the  i-th  field  of  the  following  card  contains 
the  number  of  missiles  per  launcher  for  SAM  type  i.  The 
following  card  has  a 10F8.0  format. 

The  number  of  penetrator  raid  sizes  to  be  considered  (1-10). 

The  following  card  (format  1015)  gives  the  respective  raid 
sizes. 

The  number  of  penetrator  offsets  considered  in  calculating 
the  probability  of  survival  of  a penetrator  encountering  a 
SAM  site.  Default  is  19. 

Defines  the  high  altitude  SAM  type  for  which  the  following 
cards  contain  the  number  of  SAMs  per  specified  n kilometer 
interval  (see  DEL).  The  first  card  following  gives  the 
number  of  intervals  (N)  which  are  to  be  specified  (format  15). 
The  remaining  N cards  give  the  number  of  the  interval 
(1  * 0-n  km,  2 * n-2n  km,  etc.)  and  the  corresponding  number 
of  SAMs,  one  Interval  per  card  (format  215). 

This  identifier  and  its  attendant  cards  are  similar  to  NSAMH 
and  its  associated  cards.  They  describe  the  distribution 
of  SAMs  effective  against  low  altitude  penetrators. 

The  number  of  penetrator  columns  to  be  formed  from  the  penetra- 
tors. A column  extends  parallel  to  the  direction  of  flight. 

The  number  of  penetrators  carrying  decoys.  Penetrators 
carrying  decoys  do  not  penetrate  further  after  releasing 
their  decoys. 
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PKH 


PKL 


PLCHRDY 


PSTART 


ROS 


SALSZ 


SP 


SRLD 


TERINO 


Indicates  that  the  i-th  field  of  the  following  card  gives  the 
salvo  P|^'s  for  SAM  type  1 against  the  "high  altitude"  penetra- 
tor  type  under  consideration.  The  following  card  has  the 
format  10F8.0. 

Similar  to  PKH  except  that  it  gives  the  salvo  Pj^'s  for  "low 
altitude"  penetrators. 

Indicates  that  the  1-th  field  of  the  following  card  contains 

the  probability  of  a launcher  of  a type  1 SAM  being  ready  3t 

the  time  the  raid  passes  (default  1).  Format  of  following 
card  Is  10F8.0. 

Allows  the  starting  of  a raid  at  any  of  several  distances  past 
the  FEBA.  The  value  should  correspond  to  an  integral  multiple 
of  the  interval  depth  DEL  (default  zero). 

Indicates  that  the  1-th  field  of  the  following  card  contains 

the  reliability  of  a salvo  of  an  1-th  type  SAM  (default  1.0). 

Format  of  following  card  is  10F8.0. 

Indicates  that  the  i-th  field  of  the  following  card  gives  the 
number  of  missiles  in  a salvo  of  SAM  type  1 (default  two). 

The  following  card  has  the  format  10F8.0. 

The  spacing  between  columns  of  penetrators  measured  in  kilo- 
meters. 

Indicates  that  the  1-th  field  of  the  next  card  specifies  the 
launcher  reload  time  for  SAM  type  1 (.In  secondsl.  Format  of 
following  card  is  10F8.Q. 

Indicates  that  the  1-tK  field  of  the  following  card  specifies 
for  SAM  type  1 whether  all  launchers  at  a site  are  affected 
independently  by  terrain  effects  (*  1.0)  or  whether  they  are 
affected  as  a group  (,«0)  (see  section  4).  Following  card  has 
a 10F8.0  format. 
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TITLE 


Indicates  that  the  case  title  is  given  on  the  following  card. 
The  format  of  the  following  card  is  10A8. 


TR 


VP 


The  spacing  between  rows  of  penetrators  measured  in  seconds. 
A row  extends  perpendicular  to  the  direction  of  flight. 

The  velocity  of  the  penetrators  measured  in  kilometers/sec. 


SECTION  4 


CALCULATION  OF  PENETRATOR  PROBABILITY  OF  SURVIVAL 

There  are  many  possible  models  for  determining  the  survival  probability  of  a 
given  penetrator  of  a raid  which  is  attacked  by  a SAM  site.  The  differences  are 
reducible  to  the  assumptions  on  which  the  models  are  based.  The  FIRE  model 
contains  the  following  assumptions: 

1.  A SAM  site  consists  of  one  or  more  launchers,  all  of  whose  missiles 
are  pooled  for  the  purpose  of  determining  the  number  of  salvos  available  to  the 
site. 

2.  The  SAM  site  uses  a shoot- look-shoot  firing  doctrine  against  a raid 
of  penetrators  until  all  salvos  have  been  fired  or  the  raid  leaves  SAM  coverage. 

3.  The  salvos  of  a site  are  uniformly  distributed  over  the  members  of 

the  raid. 

4.  Launchers  are  in  close  proximity,  allowing  the  site  to  be  treated  as 
a point  when  determining  the  lethal  envelope. 

5.  The  terrain  effects  (as  manifested  by  the  "probability  of  favorable 
terrain")  may  apply  independently  to  each  launcher  or  apply  to  the  site  as  a whole. 
These  two  cases  will  be  called  "terrain  independent"  and  "terrain  dependent," 
respectively.  Formulations  for  both  cases  will  be  derived. 

The  Initial  formulation  will  develop  the  technique  for  calculating  the 
expected  probability  of  survival  of  a penetrator  which  encounters  an  available 
SAM  site.  This  average  will  be  calculated  with  respect  to  the  offset  from  the 

Q 

site  at  which  the  raid  passes  the  site.  Penetrators  are  assumed  to  fly  straight 


9.  The  center  of  the  raid  is  used  to  determine  the  value  of  the  offset.  Although 
not  demonstrated  here,  the  raid  width  can  be  shown  to  be  of  little  importance  in 
the  final  result. 
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A 


and  level  at  constant  speed  throughout  the  encounter.  The  expected  value  is 
approximated  by  considering  a number  of  evenly  spaced  discrete  offsets  beginning 
with  site  overflight  (zero  offset)  and  progressing  to  the  extreme  of  the  nominal 
SAM  effective  radius,  R^.  These  offsets  will  be  designated  Oq,  0^,  Oj^, 

respectively.  There  are  actually  two  offsets  designated  by  0^  O2,  Oj^;  one 

on  either  side  of  the  site.  Thus  the  average  survival  probability  is  calculated 
relative  to  2K  + 1 offsets,  each  of  which  represents  the  survival  probability  of 
the  raid  penetrating  in  an  interval  of  offset  centered  on  some  0^  and  km 

wide.  The  probability  of  the  raid  being  in  the  interval  represented  by  offset  0^ 
is 


1 

2K  + 1 


for 


i = G 


2 

2K  + 1 


for 


K 


For  SAM  type  i , let 


N 

P 


j 


SR 


i 


PK^ 

PLR 


1 


PT 


1 


FT. 


NL, 


= number  of  penetrators  in  the  original  raid 

= probability  a given  penetrator  in  the  raid  has  survived  to 
encounter  the  site,  i.e.,  has  survived  through  the  previous 
(j  - 1)  intervals 

= salvo  reliability 

= SAM  salvo  P|^ 

3 probability  that  a given  launcher  at  the  site  is  ready 

3 probability  of  favorable  terrain  for  terrain  dependent  case 

s probability  of  favorable  terrain  for  terrain  independent  case 

s number  of  launchers  at  the  site 
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S.  , i the  number  of  salvos  that  a site  with  exactly  n launchers 
^ ,n,k 

available  (ready)  can  fire  at  a raid  with  offset 


The  probability  of  having  exactly  n of  NL  launchers  at  a site  available  is 
given  by  the  binominal  theorem  as 


/NLA  NL.  - n 

P(n)  = PLR^d  - PLR^) 


(NL.)l  n ^^'-i  • " 

ni(NL,.  --n)! 


(1) 


The  probability  of  survival  of  a given  penetrator  of  a raid  with  offset  Oj^ 
may  then  be  written 


NL. 


1 


PS{0^}  = P(n)(l  - PK.)  ’ 

n * 0 


(2) 


for  >1.  or 

NL. 

1 

PS(O^)  - Y. 

n * 0 


(3) 


for  (S^  ^ j^SR^)/NPj  < 1.  Here  NP^  is  the  number  of  penetrators  surviving  and  is 
taken  as  the  maximum  of  NPj  and  1 (it  does  not  make  sense  to  consider  less  than 
one  penetrator).  The  determination  ^ 1^  is  discussed  in  section  5.  It  is 

calculated  for  each  offset,  for  each  possible  number  of  available  launchers  n, 
n ■ 0,  ....  NL..  Since  S.  . . = 0,  equations  2 and  3 may  be  rewritten  as 

• 1 tV  » > 


PS{0. 


NL, 

(S.  „ kSR.)/NP 

} • P(0)  + ^ P(n)(l  - PK^)  J 


(4) 


n ■ 1 
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and 


^ PS{0|^} 

i 


A 

p(0)  * 2^  P(n)(l  - PK,S,_„_uSR,/NPj) 
n * 1 


(5) 


respectively. 

The  probability  of  a given  penetrator  surviving  an  encounter  over  all  possible 
offsets  is  then  simply 


K 

PS(E,,A,.I,)  . HTT  "SfV  * 2K^  S <®> 

k = 1 

where  the  notation  {E^,A^,T^}  simply  is  intended  as  a reminder  that  site  encounter 
and  availability  have  been  assumed  for  a type  i SAM  site,  as  has  favorable 
terrain. 

The  assumption  of  favorable  terrain  is  removed  in  one  of  two  ways,  depending 
upon  whether  the  terrain  dependent  or  terrain  independent  case  applies.  Proba- 
bility of  survival  for  the  terrain  dependent  case  becomes 

PS{E^,  A^}  - 1 - PT^d  - PS{E^,  A..  T^})  (7) 

For  the  independent  case,  equation  (6)  still  applies  if  PLR^  in  equation  1 is 
replaced  by  (PLR^  i.e., 

K 

PS(E,.A,1  . 

k • 1 

Introducing  the  probability  of  the  site  being  available,  PSA^,  i.e.,  deter- 
mining the  conditional  probability  that  a penetrator  survives  a site  given 
encounter,  PS{E^},  gives 
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(9) 


PS{E^}  » (1  - PSA^)  + PSA^  X PS{E^,  A^} 

Removing  the  final  condition  to  produce  the  probability  of  survival  relative 
to  a site  of  SAM  type  i,  PS^ , gives 

PS.  ^ (1  - PE^)  + PE^  X PSCE^}  (10) 

where  PE^  is  the  probability  of  encountering  the  site. 

In  the  FIRE  model  PE^  is  taken  as 

where  C * corridor  width 
w 

R » raid  width 
w 

R^  * SAM  effective  radius  as  determined  by  calculation  (R^  1 Rg»  the 
nominal  effective  radius  of  the  SAM  when  not  subjected  to  sensor 
masking,  etc.) 

For  NS^j  sites  of  SAM  type  i in  the  j-th  interval 

NS 

PS^(NS^j)  • [(1  - PE^)  + PE^  X PSS^jlE^}]  (11) 

The  product  of  the  PS^(NS^j)  over  all  SAM  types  and  over  all  traversed  inter- 
vals gives  the  overall  probability  of  survival. 

PS  ■ n n PS. (NS,.) 

J i ’ 

When  decoys  are  present  the  number  of  penetrators  becomes  NPj  + N^P^j,  where 
N^  is  the  number  of  decoys  released  and  P^  is  their  subsequent  probability  of  survival 
through  the  (J  - l)st  Interval.  Since  decoys  may  be  released  after  the  penetrators 
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carrying  the  decoys  have  been  subjected  to  attrition,  is  given  by 


where  is  the  number  of  penetrators  carrying  decoys  and  Dp  represents  the  number 
of  decoys  carried  per  decoy-carrying  penetrator. 
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SECTION  5 


EXPECTED  NUMBER  OF  SALVOS  PER  SAM  ENCOUNTER 

A.  INTRODUCTION 

Integral  to  the  calculation  of  the  survivability  of  a penetrator  in  a raid 
as  discussed  in  section  4 is  the  determination  of  the  number  of  salvos  that  may 
be  fired  at  a raid  for  a SAM  encounter  at  a given  offset.  This  section  addresses 
itself  to  that  problem.  Even  considering  the  several  simplifying  assumptions 
used  here,  the  discussion  is  somewhat  complicated  as  it  involves  both  special  and 
temporal  relationships.  It  has  been  divided  into  five  distinct  parts:  (1)  state- 
ment of  the  problem,  (2)  the  effect  of  masking  angle  on  initial  detection  of  the 
penetrator,  (3)  determination  of  the  point  of  first  possible  intercept,  (4)  deter- 
mination of  the  point  of  last  possible  intercept,  and  (5)  weighting  and  averaging 
the  results. 

B.  STATEMENT  OF  THE  PROBLEM 

Given  a SAM  site  as  described  in  section  2 and  subject  to  a constant  masking 
angle,  approximate  the  number  of  salvos  such  a site  could  fire  at  a penetrating 
raid  for  a given  offset  from  the  site. 

C.  THE  EFFECT  OF  MASKING  ANGLE  ON  INITIAL  PENETRATOR  DETECTION 

Consider  the  case  of  a direct  overflight  of  a SAM  by  a penetrator  at  altitude 
h.  The  situation  is  depicted  in  figure  3 for  identical  front  and  rear  masking 
angles.^®  The  angles  n and  a reorrttnt  the  radar  self-masking  angle  and  terrain 
masking  angle,  respectively.  The  penetrator  will  be  detected  at  a ground  range 
of  a in  front  of  the  site  and  disappear  at  a ground  range  of  a behind  the  site. 

The  SAM  must  track  and  fire  while  the  penetrator  is  visible. 


10.  Assumption  of  constant  masking  angle  is  part  of  the  statement  of  the  problem. 

11.  For  the  radar  guided  SAMs,  the  intercept  must  take  place  before  the  penetra- 
tor leaves  sensor  coverage.  For  IR  SAMs,  only  launch  must  take  place  before  the 
penetrator  leaves  coverage  as  the  IR  missile  is  autonomous  after  launch. 

26 


Figure  3.  Masking  Angle  and  Detection  Range.  The  geometry  of  the  masking  scenario 
of  the  FIRE  nwdel  is  shown  in  a plane  containing  the  center  of  the  Earth,  the  SAM 
sensor,  and  the  penetrator.  Jhe  penetratoj;;_ is  unmasked  at  an  uprange  dijtance  a 
from  the  site  and  disappears  at  a downrange.  Clearly,  relocating  the  positions 
of  the  masks  or  changing  the  sensor  height  will  affect  the  value  of  a. 


The  distance  a can  be  determined  from  the  situation  depicted  in  figure  4.  As 
in  figure  3,  h is  the  altitude  of  the  penetrator.  a represents  the  total  masking 
angle  (self-masking  plus  terrain)  and  h^.  is  the  height  of  the  sensor.  In  order  to 
allow  for  refraction  of  the  radar  signals  by  the  Earth's  atmosphere,  the  radius  of 


the  Earth  is  taken  to  be  that  of  an  earth  4/3  the  size  of  the  actual  Earth.  This 

12 

decreases  the  curvature  of  the  Earth. 

It  is  necessary  to  find  the  range,  r,  at  which  the  penetrator  emerges  from  the 
mask.  We  have 


whence  from  the  sine  formula  y can  be  determined 


12.  This  approximation  is  really  of  no  consequence  to  the  FIRE  model's  results 
because  relatively  short  ranges  are  Involved  at  low  altitude  between  the  sensor 
and  penetrator. 
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Figure  4.  Geometry  for  Computii^g  Detection  Range.  The  geometry  for  determining 
the  detection  range  r of  a penetrator  (height  h)  at  the  apex  of  angle  y by  a 
sensor  at  the  apex  of  angle  8 (hy  above  ground)  Is  shown.  The  masking  angle 
a Is  the  sum  of  self-masking  and  terrain  masking.  Rg  Is  the  radius  of  the  Earth. 
For  practical  cases,  a and  r are  very  nearly  equal. 


Again,  using  the  sine  formula 


s1n(a  * b)  , s1n(n  - Y - g - a)  , sInd 
ttg  h r r 
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(R  + h)sin6 

“ sin(a'~  &)  (detection  slant  range) 

And  finally, 

a a r cosa  (detection  ground  range) 

This  result  is  also  valid  for  flight  paths  with  nonzero  offsets.  Thus,  for 
a constant  masking  angle,  the  initial  detection  will  always  take  place  on  the 
circumference  of  a circle,  the  radius  of  that  circle  being  determined  by  the 
masking  angle  or  the  radar  range  if  masking  is  not  a factor. 

The  assumption  of  a constant  mask  angle  is  a substantial  approximation.  The 
general  problem  of  terrain  masking  is  discussed  in  appendix  A. 

D.  THE  DETERMINATION  OF  THE  POINT  OF  FIRST  POSSIBLE  INTERCEPT 

Let  the  SAM  site  be  located  at  the  origin  of  a right-handed  x,  y,  z coordinate 
system  whose  x-  and  y-axes  are  perpendicular  to  a radius  vector  from  the  center  of 
the  Earth  to  the  site.  Then  the  z-axis  is  coincident  with  the  radius  vector.  The 
lethal  envelope  of  the  SAM  site  is  modeled  as  a cylinder  of  radius  r^  whose  axis 
coincides  with  the  + z-axis.  Coaxial  with  the  lethal  envelope  is  a smaller  cyl- 
inder of  radius  r^  in  which  intercepts  may  not  take  place.  This  volume  is  referred 
to  as  the  dead  zone.  The  sensor  coverage  is  modeled  as  spherical  with  an  effective 
radius  r^.  which  is  determined  by  the  range  of  the  penetrator  at  the  first  detectable 
signal.  This  point  is  determined  either  by  emergence  from  masking  or  by  the  initial 
point  of  radar  burnthrough  when  the  signal  to  noise  ratio  first  permits  detection. 

A penetrator,  flying  straight  and  level,  enters  the  radar  and  SAM  coverage  at  alti- 
tude h as  illustrated  in  figure  5. 


13.  Designated  in  FIRE  by  the  sensor  range  (see  descriptor  ISET).  This  discussion 
describes  explicitly  the  situation  for  the  radar,  but  is  applicable  with  obvious 
modification  to  the  IR  sensor  case. 


) 

> . 


Since  the  model  is  constrained  to  considering  penetrators  flying  straight 
and  level,  the  boundaries  displayed  by  a cross  section  of  the  cylinder  (repre- 
senting the  SAM  lethal  envelope)  and  the  hemisphere  (representing  radar  range) 
at  the  ponetrator  altitude,  h,  are  of  primary  interest.  Figure  6 represents  these 
boundaries  as  seen  looking  down  at  the  z-axis  at  a cross  section  taken  perpen- 
dicular to  the  z-axis  at  altitude  h.  The  radar  range  relative  to  the  bAM  lethal 
cylinder  in  this  projection  is  given  by  ■ (r^  - h^)^. 

The  situation  normally  expected  is  the  one  presented  in  figure  6(a)  in  which 
the  horizontal  component  of  the  radar  coverage  at  altitude  h exceeds  the  lethal 
range  of  the  SAM.  However,  the  situation  shown  by  figure  6(b)  in  which  r^  exceeds 
R^  could  result  from  radar  Jamming,  small  penetrator  cross  section,  or  masking  of 
the  penetrator. 
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Figure  6. 


Cross  Section  of  Radar 
Altitude  h. 


3 


The  relationship  of  time  and  geometry  will  determine  the  interaction  of  the 
SAM  and  the  penetrators.  For  convenience,  the  penetrators  are  considered  to  fly 
parallel  to  the  y-axis.  This  makes  the  value  of  the  offset  equal  to  the  x coordi- 
nate of  the  penetrator.  It  is  also  assumed  that  the  missile  can  reach  any  point 
in  the  lethal  envelope  at  an  average  speed  v^^.  The  altitude  of  the  penetrator  is 
designated  by  h,  its  speed  by  Vp,  and  the  flight  path  offset  by  x.  The  deriva- 
tion of  first  possible  intercept  in  this  subsection  and  last  possible  intercept 
in  the  next  subsection  will  treat  a single  penetrator.  The  considerations  intro- 
duced by  multiple  penetrators  in  a raid  with  width  and  depth  will  be  discussed 
in  subsection  E of  this  section. 

Examining  figure  7 we  can  derive  an  expression  for  finding  the  point  at  which 
the  first  intercept  can  take  place  in  the  case  of  a radar  guided  SAM  (the  tech- 
nique is  identical  for  an  IR  guided  SAM).  Penetrator  detection  takes  place  at 
distance  R^.  After  a time  delay,  t^,  during  which  the  penetrator  flies  i^jVp,  the 
missile  is  fired.  The  intercept  will  take  place  at  some  point  (x,  y,  z)  so  that 
the  flight  time  of  the  missile  plus  the  delay  time  will  be 


which  must  equal  the  time  it  takes  the  penetrator  to  fly  from  the  detection  point 
to  (x,  y)  or 


Squaring  this  equation  results  in  a quadratic  in  y which  has  solutions 


Figure  7.  The  Geometry  For  Determining  The  Position  of  Earliest  Possible  Intercept. 
The  SAM  which  Is  located  at  the  origin  of  the  coordinate  system  1s  being  approached 
by  a penetrator  with  velocity  Vp  and  offset  x.  The  circles  represent  the  SAM  dead 
zone  (r^),  SAM  effective  rad1us^(r  ),  and  sensor  range  (R  ). 
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where 


tj  and  Cj 


Only  one  of  these  solutions  will  normally  satisfy  equation  (1)  as  the  other 
was  introduced  in  squaring  to  form  the  quadratic.  If  the  solutions  are  imaginary 
(both  or  neither  will  be,  a result  from  the  theory  of  equations),  they  are  of  no 
interest.  Assuming  the  solutions  are  real,  several  possibilities  exist  for  the 
solution  of  interest  (satisfying  equation  (1)),  which  will  be  called  f^i^om 

the  left  side  of  equation  (1)  is  it  obvious  that  y^^^  cannot  lie  in  front  of  the 
site  and  outside  of  radar  coverage  (y^^^  > - x^).  If  y^^^  > 0 and  places  the 

intercept  in  radar  coverage  but  outside  the  lethal  envelope,  then  the  solution  as 
it  stands  is  not  acceptable.  However,  it  indicates  the  first  possible  intercept 
will  occur  at  the  boundary  of  the  lethal  envelope.  (The  SAM  must  simply  delay 
firing  a time  greater  than  t^.)  If  y^^^  ^0  and  places  the  intercept  in  the  lethal 
envelope,  it  is  a satisfactory  solution.  This  would  automatically  similarly  be 
true  for  y^^^  < 0 were  it  not  for  the  condition  placed  on  the  SAM  that  it  may  not 
fire  at  the  penetrator  after  the  penetrator's  y value  is  less  than  zero  (the 
penetrator  has  "passed"  the  site).  Thus,  for  y^^^  < 0,  we  must  check  to  ensure 
that  at  the  time  of  SAM  launch  the  penetrator  has  not  yet  passed  the  site.^^  Other- 
wise the  solution  is  unsatisfactory  from  the  model  point  of  view. 

The  time  of  intercept  is  taken  to  be  negative  for  a positive  y,  0 for  y = 0, 
and  positive  for  a negative  y.  The  time  of  first  possible  intercept  will  be 

- -^sol/'p* 

E.  THE  DETERMINATION  OF  THE  POSITION  AND  TIME  OF  LAST  POSSIBLE  INTERCEPT 

Consider  the  schematic  representation  of  a SAM  site's  coverage  and  a penetra- 
tor flight  path  as  shown  in  figure  8.  The  point  P^  represents  the  position  of 


14.  This  restriction  Is  not  Imposed  on  IR  SAMs  which  may  fire  until  the  raid  leaves 
sensor  coverage. 
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Figure  8.  The  Geometry  for  Determining  the  Position  of  the  Last  Possible  Inter- 
cept. The  situation  displayed  Is  for  a radar  guided  SAM  as  1s  located  at 
the  point  the  penetrator  passes  the  site.  The  arguments  In  the  text  In  no  way 
depend  upon  this  choice. 


the  penetrator  at  the  Instant  the  first  salvo  may  be  fired.  This  point  Is  deter- 
mined by  penetrator  detection  range  and  the  required  time  delay  from  detection  to 
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initial  launch.  The  point  represents  the  position  of  the  penetrator  at  the 
last  instant  that  a salvo  may  be  fired.  This  point  is  determined  by  air  defense 
doctrine  or  intrinsic  geometrical  SAM  limitations  (as  opposed  to  limitations 
imposed  by  terrain,  penetrator  speed,  etc.).  The  point  lies  on  the  flight 
path  at  the  point  of  the  last  conceivable  intercept.  In  the  extreme,  P^  may  lie 
on  the  circumference  of  the  circle  defining  the  limit  of  SAM  effectiveness.  Other- 
wise, it  will  lie  closer  to  the  SAM  as  dictated  by  the  masking  angle.  For  a radar 
guided  SAM,  the  y value  of  P^^  will  have  the  value  Vb^  - x^,  where  b is  the  effec- 
tive SAM  radius  or  the  ground  range  at  which  the  penetrator  is  masked,  whichever 
is  smaller.  For  an  IR  SAM,  the  y value  of  P^  will  be  - x^,  where  b is  the 
effective  SAM  radius  or  the  y value  of  the  penetrator  at  the  instant  of  intercept, 
given  the  salvo  was  fired  at  the  time  the  penetrator  was  masked;  again,  whichever 
is  smaller. 

If  it  is  specified  that  the  missile  will  intercept  the  penetrator  at  P^ 

(regardless  of  whether  this  is  consistent  with  the  definition  of  the  SAM  capabil- 
ities), the  position  of  the  penetrator  at  SAM  launch  must  necessarily  lie 
(1)  between  P^  and  P^,  (2)  between  and  P^,  or  (3)  outside  P^.  The  significance 
of  each  of  these  locations  is  as  follows:  When  the  penetrator  lies  between  P^ 
and  P^  (Case  2),  an  intercept  can  actually  occur  at  P^.  For  the  penetrator  being 
positioned  between  P^  and  P^  (Case  1),  no  intercept  can  occur  at  P^  since  this 
would  imply  a launch  not  permitted  by  the  definition  of  the  SAM.  In  this  instance, 
the  penetrator  is  so  slow  relative  to  the  missile  that  the  latest  acceptable  inter- 
cept will  occur  when  the  missile  is  fired  when  the  penetrator  is  at  P^.  The  situa- 
tion of  Case  3 is  that  the  penetrator  is  so  fast  relative  to  the  missile  that  an 
allowable  intercept  at  P^  is  not  possible.  In  this  instance,  the  latest  acceotable 
intercept  will  occur  when  the  missile  salvo  is  fired  at  the  instant  the  penetrator 
is  at  P^.  (However,  it  is  quite  possible  that  a launch  at  this  time  will  not  result 
in  an  intercept  at  all.)  This  allows  the  missile  the  maximum  time  to  fly  to  intercept. 

In  summary,  the  above  discussion  indicates  how  to  determine  the  position  of  the 
penetrator  at  SAM  launch  which  results  in  the  latest  possible  intercept.  This  makes 
it  straightforward  to  calculate  also  the  position  of  that  intercept  and  its  time  of 
occurrence. 


15.  See,  again,  the  definition  of  P^^. 
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For  Case  1 (missile  launch  at  P^),  the  penetrator  will  fly  a distance^^ 

at  speed  v in  time  (yp^^  - -^int^^^’  missile 

will  fly 


i 


X *><  * yi„t 


Hence, 


* yfnt 


which  yields  solutions 


where 


C 


1 


2 2 
+ h 


(2) 


16.  In  the  following  discussion  the  y subscripts  have  the  following  meanings: 
int  * y coordinate  of  the  penetrator  at  intercept;  y , y , y.  are  the  y values 
r the  points  P^,  Pj,  and  A.  Pf  Pj 
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As  these  solutions  were  obtained  by  squaring  (2)  and  solvirg  the  resulting 

quadratic,  only  one  satisfies  equation  (2).  The  appropriate  y will  lie  between 

P,  and  P . The  time  of  intercept  is  -y.„*/v^. 

At  int  p 

For'' Case  2 (Intercept  at  P^) , the  time  of  Intercept  is  simply 


where  r Is  the  effective  radius  of  the  SAM  or  that  determined  by  the  masking 
e 

angle. 

Finally,  for  Case  3,  the  time  of  flight  of  the  penetrator  is 


^P^  - ^int 

which  equals  the  time  of  flight  of  the  missile 


(3) 


where 


C 


1 


Pf  2 

- and  Cj  = C-j 

P 


+ 


Again,  only  one  solution  satisfies  the  original  equation,  equation  (3).  The 
appropriate  will  lie  between  and  0,  if  a y^.^^  exists.  The  corresponding 
time  of  intercept  is  -y^^^/Vp. 

F.  THE  EFFECT  OF  MULTI-PENETRATOR  RAIDS  AND  THE  AVERAGE  NUMBER^ ^ OF  SHOTS  PER 
RAID 

Penetrator  raids  as  modeled  in  FIRE  possess  width  and  depth.  To  a high 
degree  of  approximation,  the  number  of  shots  fired  at  a raid  by  a site  is  indepen- 
dent of  the  raid  width.  However,  raid  depth  can  be  very  important  in  determining 
the  number  of  salvos  because  the  deeper  a raid,  the  longer  it  is  exposed  to  fire. 

As  a raid  passes  a SAM  site,  its  members  may  be  fired  upon  at  various  dis- 
tances from  the  site.  This  should  be  considered  in  computing  the  total  number  of 
shots.  First,  flyout  time  is  a linear  function  (in  FIRE)  of  the  site  to  penetrator 
distance  at  intercept.  Second,  where  the  raid  will  be  fired  upon  is  open  to  ques- 
tion and  can  probably  at  best  be  treated  as  a random  phenomenon.  The  basis  of  the 
method  for  handling  this  is  shown  in  figure  9.  A SAM  site  and  penetrator  flight 
path  are  depicted.  At  any  point,  F,  on  the  flight  path  where  an  intercept  may  be 
performed,  the  raid  will  take  a time  to  pass  of 

raid  depth 
penetrator  velocity 

If  it  is  assumed  that  the  front  of  the  raid  is  intercepted  at  F,  then  depending 
upon  raid  depth,  time  between  salvos,  and  flyout  time  for  the  missile,  the  raid 


17.  The  method  of  averaging  described  here  is  obviously  only  one  of  many  which 
might  have  been  used.  The  choice  of  this  one  merely  reflects  the  choice  of  the 
analyst. 
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Figure  9.  Method  of  Determining  the  Number  of  Salvos  Fired  for  a Given  Offset 
Fired 

may  be  additionally  intercepted  at  F,  0,  1,  2,  or  more  times  before  it 
passes  F (only  Integers  are  permitted).  This  number  varies  purely  as  a func- 
tion of  flyout  tin»,  as  raid  depth  and  time  between  salvos  remain  constant  as 
the  location  of  F varies. 

FIRE  makes  this  calculation  for  a nunber  of  points,  F,  evenly  spaced 
between,  and  including,  the  points  of  first  and  last  possible  Intercept.  A 
weighted  average  flyout  time  is  calculated  where  the  weighting  factors  are  the 
total  salvos  which  can  be  fired  as  the  raid  passes  a given  point.  This  average 
gives  definition  to  the  notion  that  the  raid  will  be  attacked  more  at  short 
distances  from  the  site  than  at  longer  distances.  The  number  of  shots  for  the 
given  offset  is  then  taken  simply  as 

time  of  last  intercept  - time  of  f'trst  intercept  - time  in  dead  zone 
average  flyout  time  + time  delay  between  salvos 
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where.  If  necessary,  the  exhaustion  of  missiles  on  launchers  and  subsequent 
reloading  time  are  taken  into  account  to  reduce  this  number  appropriately.  The 
calculation  is  performed  for  each  of  the  K offsets  discussed  in  section  4 for 

j • each  of  the  possible  number  of  available  salvos. 

! 
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APPENDIX  A 


THE  TERRAIN  MODELING  PROBLEM 

There  is  no  aspect  of  the  survivability  problem  for  low  altitude  penetrations 
as  important  and  as  difficult  to  treat  accurately  as  that  of  the  interaction  of 
the  terrain  with  SAMs  and  penetrators.  This  appendix  has  been  included  specif- 
ically to  discuss  this  interaction,  the  reasons  for  its  importance,  and  those 
factors  which  must  be  considered  in  any  detailed  model. 

For  penetrators  whose  altitudes  are  substantially  below  300  meters,  terrain 
effects  may  drive  survivability  results.  This  is  illustrated  in  figure  A-1  by 
results  obtained  from  the  FIRE  model.  Four  curves  giving  the  survival  probability 


Figure  A-1.  The  Effect  of  Masking  Angle  on  Survivability.  The  effect  of  masking 
angle  on  the  survivability  of  a single  penetrator  at  90  m altitude  is  depicted. 
For  each  curve,  the  constant  masking  angle  indicated  was  given  to  all  defenses. 


A-l" 


Elevation  Angle 


of  a single  penetrator  as  a function  of  penetration  distance  are  plotted  for  a 
90  tn  altitude  penetration.^®  Each  curve  represents  a different  constant  masking 
angle  about  all  SAM  sites.  The  effect  of  increasing  masking  angle  is  dramatic. 
Unfortunately  for  the  analyst,  the  ideal  situation  of  constant  masking  angle  is 
not  representative  of  the  real  world.  The  reason  for  the  pronounced  dependence 
on  masking  angle  can  be  seen  in  figure  A-2  which  shows  penetrator  elevation 
angles  for  a particular  flight  path  offset  and  several  altitudes  as  a function 
of  the  time  in  seconds  from  passing  the  site. 


Figure  A-2.  Penetrator  Elevation  Angle.  The  elevation  angles  of  various  altitude 
penetrators  with  1 kilometer  offset  from  a sensor  3 meters  above  ground  is  given 
as  a function  of  time  before  (or  after)  passing  the  site.  This  figure  and  the 
previous  figure  should  be  considered  together. 


18.  These  cases  are  typical  of  a heavy  defense  environment.  As  they  are  merely 
illustrative,  no  additional  details  of  the  scenario  are  given. 

A-2 


The  first  requirement  in  constructing  a detailed  terrain  model  is  an  under- 
standing of  the  issues  involved.  First,  it  must  be  recognized  that  the  deter- 
mination of  the  terrain  masking  around  a SAM  is  partly  a function  of  chance  and 
partly  a function  of  choice.  For  example,  a mobile  SAM  escorting  a convoy  is 
more  or  less  relegated  to  accepting  the  masking  angles  imposed  by  the  terrain  and 

vegetation  surrounding  the  road.  Such  masking  angles  are  likely  to  have  a distri- 

19 

bution  worse  than  that  of  sites  chosen  at  random.  On  the  other  hand,  premedi- 

20 

tated  site  selection  could  result  in  a much  more  favorable  distribution  of 
masking  angles  than  for  sites  selected  at  random.  Ideally  then,  a good  terrain 
model  must  include  the  option  of  selecting  different  terrain  scenarios  for  each 
type  of  SAM,  and  possibly  for  subsets  of  each  type  of  SAM. 

For  an  expected  value  model  such  as  FIRE,  the  foregoing  discussion  implies 
that  a means  must  be  devised  of  statistically  describing  the  mask  angles  about  a 
site.  A simple  example  shows,  unfortunately,  that  this  distribution  should 
ideally  also  include  information  about  distances  to  the  masking  features.  For 
example,  a distant  mask  will  hide  a close  penetrator  seen  with  the  mask  as  back- 
ground if  the  clutter  rejection  capabilities  of  the  SAM  sensor  are  not  adequate. 
Such  a statistical  distribution  must  also  take  vegetation  and  man-made  structures 
into  account.  This  can  be  seen  from  table  A-1  whicn  gives,  as  a function  of 
observer  to  mask  distance,  the  increment  in  distance  above  the  observer  necessary 
to  produce  the  specified  masking  angle.  It  can  be  seen  from  the  table  that  a 
12  meter  tree  will  constitute  approximately  a 5 degree  mask  at  100  meters 
distance  for  a 3 meter  high  observer. 


19.  Roads  tend  to  systematically  exclude  sites  such  as  hilltops  and  include 
valleys. 

20.  Sites  selected  to  provide  good  covera(^a  for  the  SAM  must  necessarily  increase 
the  visibility  of  the  SAM  to  any  penetrator.  There  appears  to  be  no  lack  of 
supporters  of  the  position  that  this  tradeoff  is  uppermost  in  the  minds  of  site 
commanders. 
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Table  A-1 

EXCESS  ABOVE  OBSERVER  HEIGHT  (m)  NEEDED  TO 
CREATE  SPECIFIED  MASKING  ANGLE 


OBSERVER  TO  MASK 
DISTANCE 
(m) 

M A 

S K I 

N G 

A N G 

L E 

r 

20 

3" 

4° 

5° 

10 

.17 

.35 

.52 

.70 

.87 

20 

.34 

.70 

1.05 

1.40 

1.75 

30 

.52 

1.05 

1.57 

2.10 

2.62 

40 

.69 

1.40 

2.10 

2.80 

3.50 

50 

.87 

1.75 

2.62 

3.50 

4.37 

75 

1.31 

2.62 

3.93 

5.24 

6.56 

100 

1.75 

3.49 

5.24 

6.99 

8.75 

200 

3.49 

6.98 

10.48 

13.99 

17.50 

500 

8.73 

17.46 

26.20 

34.96 

43.74 

1000 

17.46 

34.92 

52.41 

69.93 

87.49 

Another  major  aspect  of  the  problem  concerns  the  degradation  of  the  SAM 
capabilities  because  of  intermittent  tracking  of  the  penetrator  caused  by  masking. 
How  do  breaks  in  track  affect  shot  opportunities?  How  do  breaks  in  track  affect 
kill  probability?  What  is  the  probability  that  track  will  be  broken  at  the  time 
of  intercept?  These  are  all  potentially  important  questions  which  are  dependent 
on  the  design  of  the  SAM. 

The  FIRE  model  has  not  adequately  addressed  any  of  these  issues.  As  a result, 
the  low  altitude  survival  probabilities  may  differ  somewhat  from  those  which  would 
be  obtained  from  a more  detailed  treatment  of  the  terrain. 
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PROGRAM  LISTING 


PRCORAM  r:RE( input. output. CE8UG«CUrPUT) 


P9C0RAM  PIPE  I A A CETERMINISTIC  AIRSCRNE  PENETRATCR  SURVIVAaiL- 
ity  model  «hich  considers  Raids  CP  PENETRATCRS  against  up  TC  ’EN 
distinct  SAM  TYPES  LOCATED  WITHIN  A CORRIDOR  AS  A PUNCTICN  CP 
distance  PRCM  the  PE3A.  THE  MODEL  IS  DESCRIBED  IN  DAS-TR-7S-3. 
PIPE  A DETERM INISTC  MODEL  POR  ESTIMATING  SURVIVABILITY  THROUGH 
AREA  SURPACE  TO  AIR  MISSILE  DEPENSES.  OCTOBER  776. 


ALL  distances  in  PROGRAM  PI  RE  ARE  EXPRESSED  XILOMETERS. 
penetration  DISTANCE  IS  CONSIDERED  PLUS  3EHI\  THE  PEBa . MINUS 
IN  PRCNT  OP  the  PEBA, 

COMMON  /A/  NPHI . ITYPE.OT.RS.VM.HM.VP.RDE?.NP.DTP 
CCMMON  /a/  DZONE.RMAX 

COMMON  /C/  ILCW. ALPHA! 10) .HR. ISET( 10) . ISETL( 10) .ASPECT! IC: 

COMMON  /D/  DETCRA.PI 

COMMON  /E/  PKX.PTX.NLCH.XPPX.PSSEAIU.OEPTH.AVSHT 
COMMON  /P/  NO  I ST 

COMMON  /G/  TO( 10) .SRI  10) .SV( 10) .NS( 10) .2Z( 10) .HRS! 10) .RR( 10) .PD( iO 
1 ) .PT( 10) .PKSSI 10) .TDPI 10) 

COMMON  /J/  ENVLl 1 10) .PEN! 10) .XSITES! 10) 

COMMON  /K/  NMLI 10) .PDLI 10) .PTli 10) .TDL! 10) .SRLI 10) .SVL( 10) .DZLt 10) 
1 .RRLI 10) .HRSLI 10) .PKSLI 10) . ITITLEI 1 1 ) .TDLP! 10) .NLCHRLl IC) . TERINDl 1 
SO) 

COMMON  /PS/  NPENT.NPTI 10) .SALSZI 10) .NLCHRl 10) .PLCHRCYI 10) .RCSl 10) 
COMMON  / P3/  PENC . TR . SP . CCRUO TH ! 2 ) . DEL 

COMMON  /PH/  PSTART. ICHAPP.CHAPPS.CHAPPE, lECM.ECMS.ECME. IDECOY.DSTR 
I .DENC.CHEPPi 10) .ECMEPPI 10) 

COMMON  /PHH/  IPINT.DFIIO) .DECPPENIIO) .PENCDEC 
COMMON  /NHLN/  NSAM(66. lO.S) 

DIMENSION  ENI 10) . ITP(IO).  PSilO) 

DIMENSION  SPPXdO).  DECSTRTllO) 

DIMENSION  XSITS(IO) 

CALCULATE  DEGREES  TO  RADIANS  CONVERSION  FACTOR 
PI-H.*ATANI 1 . 1 
DETCRA-PI/ieO. 
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BESI  AVAILABLE  (DPlf 


ooo  o oo  oo  o o o o ooo  o o oor» 


aeciN  LOOP  ON  cases 

2 CONTINUE 
CALL  INPUT  I 

SELECT  HIGH  OR  LOW  DEFENSES 
I alt -2- 1 LOW 
PRINT  CASE  title 
PRINT  13.  (ITITLE(l). 1-1.8) 

BEGIN  LOOP  ON  penetration  RAID  SIZE 

00  m M-l.NPENT 

SET  THE  total  number  OF  PENETRATORS  IN  THE  RAID 
NP-NPT(M) 

CALCULATE  the  NUMBER  OF  ROWS  OF  PENETRATCRS 
NROW-NP/ PENC -0 . 99999 

CALCULATE  THE  DEPTH  AND  WIDTH  OF  THE  RAID  IN  KM 
DEPTH- ( NRCW-1 ) -TR* VP 
W10TH-(PENC-1 . )»SP 

PRINT  RAID  information 
PRINT  IT  NP  width  DEPTH 

initialize 'running  PENETRATOR  and  DECOY  survival  FRC8AaiLI"IES 
FOR  CURRENT  INTERVAL  AND  ALL  INTERVALS  TO  DATE 
SURVP- 1 . 0 
SURVPT-1 ,0 
SURVD-l .0 
SURVOT-1 .0 

initialize  the  number  of  decoys  surviving  to  release  (ZERO  value 

USED  AS  TEST) 

00  3 I-I.IO 
CECSTRT(I)-0 

3 CONTINUE 
INITIALIZATION 

00  H 1-1,10 

XSITS(n-0. 

PS( I )-l .0 

H continue 

POIST-PSTART 
INNl-0 

BEGIN  LOOP  ON  PENETATION  DISTANCE  INTERVAL 

00  13  J-t.NOIST 
C CALCULATE  PRESENT  PENETRATION  DISTANCE 

POIST-POIST-OEL 

C CALCULATE  PRESENT  PENETRATION  INTERVAL 

lOIST-POIST/OEL 

C initialize  IX  TO  INDEX  SAM  RESULTS 

IX-0 


BEST  AVAIIABIE  COPY 


oo  o 


BEGIN  LOOP  CN  SAMS  TO  DETERMINE  PRC0ABILITY  OF  SURVIVAL  VALUES 
00  12  I-l.lO 

C DETERMINE  IF  SAM  TYPE  I IS  TO  BE  CONSIDERED  FOR  THIS  CASE 

IF  (DF(  I ) .EQ.O.OR.  ( ILOW.f.Q.C.ANO.  ISETl  I ) .EQ.Q)  .OR.  i ILCW.EQ.  I .AND.  I 
ISETLU  ) .EQ.O)  ) GO  TO  12 
C SAM  TYPE  I IS  TO  BE  CONSIDERED 

IX-IX*I 

C CALCULATE  THE  NUMBER  CF  SITES  OF  SAM  TYPE  I J THE  J-TH  BAND 

XSITES( IX)-NSAM( lOIST^l . I . IALT)-NSAM( IDIST. I . I.  ) 

C STORE  SAM  TYPE  FOR  LATER  PRINTING 

ITP( IX)-I 

C SELECT  IF  HIGH  OR  LOW  ALTITUDE  SAM  PARAMETERS  ARE  TO  BE  USED 

IF  ( ILCW.EQ. n GO  TO  5 

C CALCULATE  SAM  AVAILABILITY  AND  SET  TERRAIN  DEGRACATION  AND  ?K 

C FOR  HIGH  altitude  SAM  SALVOS 

?OX-PD( I ) ‘OFI I ) 

PTX-PT(  I ) 

NLCH-NLCHR ( I ) 

PKX»PKSS(  I ) 

GO  TO  6 

C CALCULATE  SAM  AVAILABILITY  AND  SET  TERRAIN  DEGFAQAT ICN  AND  FK 

C FOR  LOW  altitude  SAM  SALVOS 

5 CONTINUE 

POX-PCL( 1 ) •CF( I ) 

PTX»PTL( I ) 

NLCH-NLCHRL  t I ) 

PKX-PKSL ( I ) 
a CONTINUE 

C CHECX  FOR  PRESENCE  CF  CHAFF  IN  THIS  INTERVAL 

IF  ( ICHAFF.NE. I .OR.PDIST.GT.CHAFFE.CR.PDIST-DEL.LT.CHAFFS)  CO  TO  7 
C CHAFF  IS  PRESENT  IN  INTERVAL.  SET  CHAFF  DEGRACaTICN 

XM«CHEFF{ I ) 

00  TO  9 

C CHAFF  IS  NOT  PRESENT  IN  INTERVAL 

7 CONTINUE 
XM-1 .0 

C CHECK  FOR  PRESENCE  OF  ECM  IN  THIS  INTERVAL 

9 CONTINUE 

IF  ( lECM.NE.  I ,CR.PO!ST.GT.EC,M£.OR.PO!ST-OEL.LT.ECMS)  OC  TO  3 
C ECM  PRESENT  IN  INTERVAL.  REVISE  SAM  SALVO  PK 

PKSS 1 -ECMEFF  (I)Ml.O-(l.O  -PKX 1 • • ( 1 . 0 / SALSZ  U ) 1 ) 

PKX-I .-( I ,-PKSSl )»»SALS2( I ) 

9 CONTINUE 

C CHECK  FOR  PRESENCE  OF  DECOYS  IN  THIS  INTERVAL 

IF  ( IDECOY.NE. I .CR.POIST.GT.DENO.OR.PCIST-CEL.LT.OSTR)  CO  TO  10 
C DECOYS  ARE  PRESENT  IN  INTERVAL 

C CALCULATE  INITIAL  NUMBER  OF  DECOYS  IF  FIRST  INTERVAL  WITH  DECOYS 

IF  (DECSTRTt I) ,EQ.O»  OECSTRT ( I ) -CECPPENt I ) •SURVPT*PENCDEC 
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BKI  AVAlLiBLE  COPY 


oo 


C CALCULATE  NUMBER  OF  PENETRATORS 

XPPX-SURVPTmnP-PENCOEC)^SURVOT»OECSTRT(  1 ) 

IF  (XPPX.lt. I .0)  XPPX-l.O 
IF  (PENSUR.lt. I .0)  PENSUR-l.O 
GO  TO  11 

C DECOYS  ARE  NOT  PRESENT  IN  INTERVAL.  CALCULATE  NUMBER  OF 

C PENETRATORS 

10  CONTINUE 
XPPX»SURVPT*NP 

IF  ( IDECOY.EQ. I .ANO.POIST.GT.DENO)  XPPX-SURVPT*  ^-PENCOEC) 

IF  (XPPX.LT.  1 .0)  XPPX-1.0 

11  CONTINUE 

C SAVE  total  PENETRATORS  SEEN  BY  SAM  TYPE  I (INCLUOING  OECCYS) 

SPPX( IX)-XPPX 

IF  (J.EQ. I .OR.XSITES( IX) .NE.O)  CALL  NSHOTS  (I) 

IF  (J.EQ. I ) EN(IX)-AVSHT 

C calculate  PROBABILITY  A PENETRATCR  SURVIVES  A SINGLE  SITE  OF 

C type  I IN  THE  JTH  BANO.  GIVEN  SITE  ENCOUNTER 

PSSEIJ-l .0-POX*POX*PSSEAIJ 

CALCULATE  probability  OF  ENCOUNTER  WITH  TYPE  I DEFENSES  IN  THE 
JTH  BANO 

PEN( IX)-AMINl (B.O* (ENVLl ( I ) fWlDTH) /CCRWOTH( lALT) . 1 .C) ‘XM 
CALCULATE  THE  PROBABILITY  OF  SURVIVING  ONE  SITE 
PSl IJ-I .0-PEN( IX)^PEN( IX)»PSSEIJ 
C CALCULATE  the  PROBABILITY  OF  SURVIVING  ALL  SITES  CF  SAM  TYPE  I 

C IN  THE  interval 

PSNIwi-PSl  1J«*XSITES(  IX) 

C calculate  survival  probability  of  PENETRATCR  RELATIVE  TO  I-TH 

C type  SAM 

PS( IX)-PS{ IX)»PSNIJ 

C CALCULATE  PENETRATCR  SURVIVAL  PROBABILITY  TO  PRESENT 

SURVP«SURVP»PSNI J 

c calculate  decoy  survival  probability  to  present  if  applicable 

IF  ( IDECOY.EQ. I .ANO.POIST.LE.CENO.ANO.POIST-OEL.GT.DSTR.ANO.CECSTR 
it:!). NE.O)  SURV0-SURV0*PSNIJ 
XS1TS(  !X)«XSITS(  IX)i-XSITES(  IX) 

la  continue 
survpt-survp 
SURVOT-SURVO 

C calculate  the  number  of  PENETRATORS  SURVIVING 
pensur-survpt»np 

IF  ( IDECOY . EQ . I . ANO . PO I ST-OEL . GE . 0S7R ) PENSUR-SURVPT* ( NP-PENCDEC ) 

C increment  print  interval  COUNTER 

INNI«INNI>1 

C determine  if  OUTPUT  DUE  THIS  INTERVAL 

IF  (INNl.NE.lPINT)  00  TO  13 
C RESET  PRINT  INTERVAL  COUNTER 

lNNl-0 

C OUTPUT  SURVIVAL  PROBABILITIES 

PRINT  16,  POIST.PENSUfl.SURVP. ( ITP( I ) .EN( I ) .PEN( I ) .XSITSI 1 ) .SPPX( I ) 


BtSl  AVAILABLE  COPY  b-a 


ooo 


I 


p 

1 .PS( I) . I-l . IX) 

13  CCNTINtt 
IH  CONTINUE 

C CHECK  FOR  AOOITIONAL  CASES 

GO  TO  2 


15  FORMAT  ( 1HI///1OX.0AIO) 

16  FORMAT  (//19H  PENETRATION  0EPTH-.F5. I .5X.32HNUM'  < OF  PENETRATORS 
1SURVIVINC-.F5. I .5X.H2HPR08ABILITY  OF  GIVEN  PENl  ATOR  SURVIVING-. F 
2^.2/ IX.  HHTYPE . 6X . 7HSAL  VOS/  . 7X . 5H  PROS . I OX . SHTO*.  . / . 5X . 1 1 HPENETRA  TO 
3RS . 5X . 8HSURV I VAL/  1 X . ‘♦HS I TE . 7X . 4HS I TE . 9X . 9HENC0UNTER . 7X , 5HS I TES . 8X . 
‘4‘+HSEEN,7X.1  lHPRC8ABILlTY/(  I3.2Fm.2.r  12.  I .2Fm.2) ) 

17  FORMAT  (//23H  NUMBER  OF  PENETRATORS- . 15/ . 12H  RAID  WIDTH-.F6.2/ :2H 
IRAIO  0EPTH-.FS.2// ) 

END 


I 
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SUBROUTINE  NSHOTS  (D 


subroutine  NSHOTS  IS  CALLED  FROM  PROGRAM  FIRE  ANO  IN  TURN  CALLS 
SUBROUTINE  SHOTS.  THE  PURPOSE  OF  NSHOTS  IS  TO  1)  ESTABLISH  THE 
values  of  PARAMETERS  USED  IN  SUBROUTINE  SHOTS  (WHETHER  FOR  THE 
HIGH  OR  LOU  altitude  SAM  SCENARIO!  ANO  2)  TO  A 'PAGE  THE  RESULTS 
RETURNED  FROM  SHOTS.  SHOTS  IS  CALLED  FOR  A SERI  S OF  SINGLE 
COLUMN  FLIGHT  PATH  OFFSETS  (RAID  DEPTH  ONLY  I'  :ONSIOEREO)  ANO 
FOR  EACH  OFFSET  RETURNS  THE  NUMBER  OF  SAM  SAl  S WHICH  MAY  BE 
FIRED  AT  THE  RAID.  THESE  NUMBERS  ARE  AVERAGEL  0 DETERMINE  THE 
EXPECTED  number  OF  SHOTS  PER  RAID  ENCOUNTER  WITH  THE  SITE.  THE 
EFFECTIVE  RADIUS  OF  THE  SITE  IS  ALSO  A PRODUCT  OF  THE  RESULTS 
obtained  from  SHOTS. 

COMMON  /A/  NPHI . ITYPE.DT.RS.VM.HM.VP.RDEP.NP.OTF 
COMMON  /8/  OZONE. RMAX 

COMMON  /C/  ILOW.ALPHAI 10) .HR. ISET( 10) . ISETLI 101 .ASPECT( 10) 

COMMON  /Cl/  SRLOdO) 

COMMON  /C2/  TRLD 
COMMON  /O/  OETORA.PI 

COMMON  /E/  PKX.PTX.NLCH.XPPX.PSSEAIJ.OEPTH.AVSHT 
COMMON  /G/  TOC  10) .SR( 10) .SV( 10! .NS( 10) .02( 10) .HRS( 10) .RR( 10) .PO( 10 
I ) .PTC iO) .PKSS( 10) .TDF( 10) 

COMMON  /I/  NX 

COMMON  /J/  ENVLl ( 10) .PEN( 10) .XSITESC 10) 

COMMON  /K/  NML(101.POL(10).PTL(10).TDL(10).SRL(10).SVL(10).OZL(10) 

1 .RRL( 10) .HRSL( 10) .PKSL( 10) . ITITLEC 1 1 ) .TOLFC 10) .NLCHfiLC 10) .TERINO: I 
20) 

COMMON  /P2/  NPENT.NPTI 10) .SALSZC 10) .NLCHR( 10) .PLCHROY( 10) .ROS( 10) 

DIMENSION  PS0K(25).  PMLRS(25) . NXS(25) 

INITIALIZE  COUNTER  OF  SALVOS  FIRED  FOR  ALL  OFFSETS.  SAM  TYPE. 
RAID  depth.  ANO  SAM  RELOAD  TIME 
XSHOT-0. 

I type- I 

rcep-depth 

TRLO-SRLD(I) 

PLR-PLCHROY( I ) 

IF  (TERINO(l).EQ.l)  PLR-PLR«PTX 
SELECT  SAM  PARAMETERS  FOR  HIGH  OR  LOW  ALTITUDE  SCENARIO 
IF  (ILOW.NE.l)  GO  TO  I 

parameters  are  for  low  altitude  scenario.  DEFINITIONS  FOR 
variables  ON  RIGHT  SIDE  OF  REPLACEMENT  STATEMENT  MAY  BE  INFERRED 
FROM  SUBROUTINE  INPUT  I ANO  DAS-TR-76-9 
RS-SRL(I) 

VM-SVL ( I ) 

DZONE-OZL(I) 


BEST  'AVAIUBIE  COPY 


ooo  o o o o oo  ooo 


Hfi-HRSL ( / ) 

0T«TDL( I ) 

0TF-T0LF( I ) 

RTEMP-RRL ( 1 ) 

GO  TO  a 

parameters  are  for  high  altitude  scenario,  definitions  for 
variables  on  right  side  of  replacement  statement  may  be  inferred 
FROM  subroutine  INPUT  I AND  OAS-TR-76-9 
I CONTINUE 
RS«SR( I ) 

VM-SV ( 1 ) 

OZONE-OZ( I) 

HR-HRS ( I I 
DT-TO( I ) 

OTF-TOF(I) 

RTEMP-RR(  I ) 
a CONTINUE 

CALCULATE  PROBABILITY  OF  EXACTLY  M LAUNCHERS  READY . INCLUDE 
terrain  if  INDEPENDENT  FROM  LAUNCHER  TO  LAUNCHER 
NLCHP-NLCH<»1 
TFAC I -NLCHP 
00  3 N-l. NLCHP 
M-N-1 
TFACa-M^l 
TFAC3-TFACI-M 

rMLrtSlN) •( iFACT(NLCH)/ ( IFACT(Mi • IFACT (NLCH-M) ) ) • (PLR*l .QE-IOO) ••M* 

I « I .0-PLR^I  .OE-IOOM*  (NLCH-M) 

CALCULATE  NUMBER  OF  AVAILABLE  SALVOS 
NXS ( N ) -M^NML ( I ) / SALSZ ( I ) 

3 CONTINUE 

INITIALIZE  OFFSET  AND  COUNTER  OF  OFFSETS  RESULTING  IN  INTERCEPTS 

z*o. 

MPOINT-0 

initialize  offset  survival  probabilities 
00  K-i.as 

PSOK(K)-0 
H CONTINUE 

determine  separation  of  ADJACENT  FLIGHT  PATH  OFFSETS 
DELTA-RS/ (NPHI-l ) 

BEGIN  LOOP  ON  PENETRATOR  OFFSET 

ZTEST-O 
00  6 J-l.NPHI 
C CALCULATE  OFFSET  TO  BE  CONSIDERED 

Z-(J-l )»DELTA 

C BEGIN  LOOP  FOR  EXACTLY  N SITES 

XSHOTRmO 
NLCHP-NLCH^ I 
00  9 N- I, NLCHP 
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PMLR«PMLRS(N) 

NX«NXS(N) 

CALL  SHOTS  (XSHOTl .Z. J.RTEMP) 

XSHOTR-XSHOTR+XSHOT 1 *PMLR 
TEST»XSHOTl«ROS( I )/XPPX 

ir  (TEST.GE. l .0)  PSOK(J)-PMLR*(l .0-PKX)«*TEST+PSOK(Ji 
IF  (TEST.LT.  1 .0)  PSOK  ( J) -PMLR*  ( I . 0-PKX*TEST  ) *PSC'<  ( J) 

5 CONTINUE 

IF  (Z.EQ.O.ANO.XSHOTR.OT.O)  ZTEST-1.0 
IF  (XSHOTR.GT.O)  MPOINT-MPOINT* 1 
XSHOT-XSHOTi-XSHOTR 

6 CONTINUE 

ALL  FLIGHT  PATH  OFFSETS  UHICH  COULD  PERMIT  INTERCEPT  HAVE  BEEN 
EXAM  I NED 

CALCULATE  THE  INTEGRATED  SURVIVAL  PROBABILITY 
PSSEAIJ»( 1 .0/ (2.0*NPHI*n ) ♦PSCK(  I ) 

BO  7 J-a.NPHI 

PSSEAIJ-(2.0/(2.0*NPHI*I .0) )«PSOX(J)*PSS£AIJ 

7 continue 

IF  (TERINOt I ) .EQ.O)  PSSEAI J«1 .0-( 1 .0-PSSEAI J) ‘PTX 
C CALCULATE  HALF-WIDTH  OF  LETHAL  ENVELOPE 

ENVL I ( I) -0 

IF  (MPOINT.NE.O.AND.ZTEST.EQ.O)  ENVLl ( I )-MP0INT»0ELTA 
IF  (MPOINT.NE.O.ANO.ZTEST.EQ. 1 .OJ  ENVLI ( I )-(MPOINT-l .0^0.5)*OELTA 
C RESET  NOMINAL  SAM  RANGE  IF  MODIFIED  IN  SHOTS 

C CALCULATE  AVERAGE  NUMBER  OF  SALVOS  FIRED  GIVEN  AN  ENCOUNTER 

IF  iMPOInT.EQ.O)  AVSHT-0 

IF  (MPOINT.NE.O.ANO.ZTEST.EQ. 1 ) AVSHT-XSHOT/ (MPOINT-l .0>0.5) 

IF  (MPOINT.NE.O.AND.ZTEST.EQ.O)  AVSHT-XSHOT/MPOINT 
RETURN 


ENO 
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subroutine  shots  (XSHOTX.X.K.RSAVE) 


c 

C SUBROUTINE  SHOTS  IS  CALLED  BY  SUBROUTINE  NSHOTS.  IT  CALCULATES 

C FOR  A SINGLE  COLUMN  FLIGHT  PATH  OFFSET  (RAID  DEPTH  ONLY  IS 

C CONSIDERED)  THE  NUMBER  OF  SALVOS  WHICH  THE  SAM  TYPE  UNDER 

C CONSIDERATION  CAN  FIRE  AT  THE  RAID.  IN  THE  PROCESS  IT  IS 

C NECESSARY  TO  determine  THE  LOCATION  OF  THE  FIRST  AND  LAST 

C POSSIBLE  INTERCEPTS.  IN  COMPUTING  THE  NUMBER  0.  SALVOS.  AN 

C AVERAGE  IS  taken  OVER  ALL  POSSIBLE  INTERCEPT  PL  NTS  WITH  SHORT 

c flyout  times  weighted  more  heavily  than  longf"  jnes.  launcher 

C RELOAD  time  IS  CONSIDERED. 

C 

C THE  Y value  of  the  PENETRATOR  IS  POSITIVE  BEFORE  PASSING  THE  SAM 

C AND  NEGATIVE  AFTER.  THE  ALGEBRAIC  SIGN  OF  TIME  IS  OPPOSITE  THAT 

C OF  THE  Y PENETRATOR  COORDINATE.  ASPECT ( I ) ,NE . 0 INDICATES  AN  IR 

C GUIDED  SAM. 

C 

COMMON  /A/  NPHI . ITYPE.DT.RS.VM.HM.VP.ROEP.NP.OTF 
COMMON  /B/  OZONE. RMAX 

COMMON  /C/  ILOW.ALPHAI 10) .HR. ISET( 10) . ISETLI 10) .ASPECT! 10) 

COMMON  /C2/  TRLD 
COMMON  /O/  OETORA.PI 
COMMON  / 1 / NX 
C 

DIMENSION  Y(2) . YFST(50).  TFSTtSO).  YLSTtSO).  TLST(50) 

C 

DATA  (ASPECT! I ).I«l.9)/8‘0.0.785H/ 

C 

C TEST  FOR  ZERO  OFFSET 

IF  (X.NE.O)  GO  TO  2 

C FOLLOWING  INITIALIZATION  DONE  ONLY  FOR  ZERO  OFFSET 

RDZ -OZONE 
TO-OT 
TOF-OTF 
H-HM 
RO-ROEP 
RE-RS 
RR-RSAVE 

C CALCULATE  MAXIMUM  LINE-OF-SIGHT  RANGE  OF  SENSOR  FOR  PENETRATOR 

C WHEN  MASKING  CONSIDERED.  USE  ^/3  EARTH  FOR  RADAR  GUIDED  SAMS 

REARTH-8H99 . 0 

IF  (ASPECT! I TYPE) .NE.O)  REARTH-837 1 . 0 

beta- AS  I N ( REARTH/  ( REARTHi-Hfl ) ) -ALPHA ! I TYPE ) ‘OETORA 

OAMMA-AStN(SIN(BETA)MREARTH-HR)/  (REARTH-H)  ) 

THE  TA -P I - ! BE  T A -GAMMA ) 

RMAX-S I N ( theta ) • ( REARTH-H ) / S I N ( BETA ) 

C DOES  RADAR  MASKING  REDUCE  SENSOR  RANGE 

IF  (RMAX.LT.RR)  RR-RMAX 

c initialize  positions  and  times  of  intercept  POINTS  OF  interest 

C -ICES  IS  A FLAG 
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00  i 

YFST(J)  — I0E6 

TFST(j)— loee 

YLST(J)  — 10E6 

tlst(J)— loee 

1 CONTINUE 

C TEMPORARILY  REDUCE  NOMINAL  EFFECTIVE  RANGE  OF  SAM  IF  NECESSARY 

IF  (RE.OT.RR)  RE-RR 

2 CONTINUE 
XSHOTX>0 

C IF  OFFSET  GREATER  THAN  EFFECTIVE  RADIUS.  RETl 

IF  (X.OT.RE)  RETURN 

C IS  FLIGHT  PATH  THROUGH  DEAD  ZONE 

IOZ-0 

IF  (X.LT.ROZ)  IDZ-l 

C DEFINE  POSITION  OF  PENETRATOR  AT  TIME  OF  LATEST  ALLOWED  MISSILE 

C LAUNCH 

YTLST-0 

IF  (ASPECT!  ITYPEJ.NE.O)  YTLST— S(a9T(RE**a-X»*2) 


SOLVE  FOR  Y VALUE  OF  FIRST  INTERCEPT 


ooo  ooooooooo 


r 


ir  (Y(J) .GT.SQRT(YRR»«2-X»*a) . AND. Y( j) .GE.O)  GO  TO  H 

C ROOT  lies  within  RADAR  COVERAGE  OR  IS  NEGATIVE 

IF  IY(J) .lT.SQRT(RE»»2-X«*2J .0R.Y( J) .LE.O)  GO  TO  3 

C FIRST  intercept  AT  BOUNDARY  OF  LETHAL  ENVELOPE 

YFST(K)-SQRT(RE**2-X»«2) 

TFST(K)— YFST(X)/VP 
TFLYM"SQRT(X»»2*YFST(K)  ••2«'H»«2)  /VM 
YPENL-YFST ( K ) ♦TFLYM* VP 
GO  TO  5 
3 CONTINUE 

C DOES  ROOT  LIE  IN  LETHAL  ENVELOPE 

IF  (AeS(Y(UJ J .GT.SQRT(R£>«2-X‘«2» ) CO  TO  ^ 

C ROOT  lies  in  lethal  ENVELOPE 

C CHECK  TO  SEE  THAT  SHOT  FIRED  BEFORE  PENETRATOR  PASSED  YTLST 

TFLYM«SQRT(X»*2*Y( J) ••2^H«»2) /VM 
YPENL-Y ( J ) ♦TFLYM»VP 
IF  (YPENL.LT. YTLST)  GO  TO  H 

C missile  fired  before  PENETRATOR  PASSED  YTLST.  SOLUTION  OK 

YFSTJK)«Y(J) 

TFST(K)«-YFST(K)/VP 

GO  TO  5 

CONTINUE 

C THERE  IS  NO  FIRST  INTERCEPT 

GO  TO  12 

C SET  Y-VALUE  of  PENETRATOR  AT  TIME  MISSILE  LAUNCHED  FCR  FIRST 

C INTERCEPT 

5 continue 

ytfst-ypenl 


FIND  POSIT ICN  AND  TIME  OF  LAST  INTERCEPT 


CALCULATE  POSITION  OF  PENETRATOR  AT  INSTANT  MISSILE  IS  FIRED 
PRODUCING  AN  INTERCEPT  AS  PENETRATOR  LEAVES  COVERAGE 
TTEST-SQRT(RE**2*H**2)/VM 
YTEST— SORT(RE**2-X**2)*TTEST*VP 
IF  IYTEST.LT.YTLST)  go  TO  8 
IF  (YTEST.LE.YTFST)  go  TO  11 

PENETRATOR  IS  TOO  FAST  FOR  LEGITIMATE  INTERCEPT  AT  COVERAGE 
AT^YTFST^^^*^  intercept  OCCURS  FOR  MISSILE  LAUNCH  WITH  PENETRATO 

CI-YTFST/VP 

A«( I .0/VP**2-t .0/VM«*2) 

B— 2.0»Cl/VP 

C«Cl**2-(X**2^H*»2)/VM**2 
TEMP-8* •2-W.0*A*C 
Y(I)-(-8^SQRT(TCMP))/(2.0*A) 
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Y(a)-(-0-SQRT(TEMP) )/ (2.0*A) 

00  6 J-1.2 
L-J 

YT-SQRTtX^’a^YCJJ’^a^H^’aJ/VM-YTFST/VP^rcjl/VP 
IF  (A8S(YT) .LT.O.OOI)  GO  TO  7 

6 CONTINUE 
STOP  5 

7 CONTINUE 
YLST(K)-Y(L) 

TLST(K)«-YLST{K)/VP 
GO  TO  12 

8 CONTINUE 

c penetrator  is  too  slow  for  legitimate  intercept  at  coverage 

C LIMIT,  latest  intercept  OCCURS  FCR  MISSILE  LAUNCH  WITH  PENETRATO 

C AT  YTLST 

CI-YTLST/VP 

A»( I .0/VP**2-I .0/VM«*2) 

8»-2.0*CI/VP 

C«C1»»2-(X*«2>H*»2)/VM**2 
TEMP-8* •2-H.0»A«C 
Y< I)-(-8*SQRT(TEMP) )/ 12.0*A) 

Y(2)-(-8-SQRT(TEMP) )/ (2,0*A) 

00  9 v'-l,2 
L-J 

YT-SQRTJX**2*Yrj)**2>H**2)/VM+Y(J)/VP-YTLST/VP 
IF  (A0S(YT) .LT.O.OOI)  GO  TO  10 

9 CONTINUE 
STOP  10 

10  CONTINUE 
YLST(K)-Y(L) 

TLST(K)— YLST(K)/VP 
GO  TO  12 

11  CONTINUE 

C LAST  INTERCEPT  TAXES  PLACE  AT  COVERAGE  LIMIT 

YLST(X)— SCRT(RE**2-X**2) 

TLST(K)— YLST(K)/VP 

c 

C 

c 

C FINO  weighting  factors  for  specific  values  of  T.  TFST(K) .LE.T.LE 

C .TLST5KI 

c 

C 

c 

12  CONTINUE 

C TEST  FOR  VALIO  INTERCEPTS 

IF  (TFST(K) .Ne.-1QE8)  GO  TO  13 
C NO  INTERCEPTS  P0SSI8LE 

XSHOTX-0 
RETURN 
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c MISSILES  WILL  BE  EXHAUSTED.  IS  THERE  ENOUGH  TIME  TO  RELOAD 

IF  (NXMTFLYA>TD)>TRLD.LE.  (TLST(K)-TrST(K)  ) ) GO  TO  19 
0 RAID  WILL  PASS  BEFORE  LAUNCHERS  CAN  BE  RELOADED 

XSHOTX-NX 
RETURN 

19  CONTINUE 

C CALCULATE  NUMBER  OF  FIRE-RELOAD  CYCLES  PERMITTED  BY  AVAILABLE 

C TIME  AND  SUBSEQUENT  NUMBER  OF  SALVOS  WHICH  CAN  FIRED 

NCYCLE-XCYCLE-(TLST(I<)-TFST(K))/(NXMTFLYA+TD)>TF>  D) 
FCYCLE-XCYCLE-NCYCLE 

IF  (FCYCLE^lNX^ITFLYA-t-TDl-fTRLD)  .GT.NXMTFLYAfTI  GO  TO  20 
XSH0TX-NCYCLE*NX>FCYCLE* ( (NX* (TFLYAfTD) ♦TRLD) / i .‘TFLYA^TD) ) ‘NX 
RETURN 

20  continue 

xshotx-(ncycle*i)*nx 

RETURN 

C 

END 


1 

! 

1 
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FUNCTION  IFACT  (I) 

FUNCTION  IFACT  CALCULATES  I FACTORIAL 
IFACT-l 

IF  (I.EQ.O)  RETURN 
00  1 J-l.I 
IFACT-IFACT»J 
1 CONTINUE 
RETURN 
END 
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subroutine  input  I 


subroutine  input  I handles  all  input  for  program  fire.  INPUTl  IS 
called  whenever  fire  requires  NEW  DATA  FOR  ADDITIONAL  CASES.  THE 
BASIS  FOR  INPUT  IS  THE  ALPHANUMERIC  IDENTIFIER  (SEE  I COM  ARRAY) 
WHICH  IDENTIFIES  THE  INPUT  VALUE  FOLLOWING  IT  ON  THE  CARO  AND/OR 
INDICATES  THE  DISPOSITION  OF  THE  INPUTS  ON  SUCCESSIVE  CAROS.  IN 
GOING  FROM  CASE  TO  CASE  ONLY  INPUTS  WHICH  CHAN  " MUST  BE  GIVEN 
(UNCHANGED  INPUTS  ON  THE  SAME  CAROS  MUST  BE  REF  ATED) . THE 
ORDER  OF  inputs  IS  IMMATERIAL  EXCEPT  WITHIN  f'  i OF  CAROS 
ASSOCIATED  WITH  AN  IDENTIFIER.  UNRECOGNIZED  NTIFIERS  PREVENT 
RETURN  TO  FIRE  FROM  INPUTl.  INPUTS  DISCUSSED  . OAS-TR-76-g 

COMMON  /A/  NPHI . ITYPE.DT.RS.VM.HM.VP.RDEP.NP.OTF 
COMMON  /B/  OZONE. RMAX 

COMMON  /C/  ILOW.ALPHA(  10)  .HR.  IS£T(  10)  . ISETL(  !(])  .ASPECT(  10) 

COMMON  /Cl/  SRLDdO) 

COMMON  /O/  DETORA.PI 
COMMON  /F/  NO  I ST 

COMMON  /G/  TO( 10) .SR( 10) .SV( 10) .NS( 10) .DZ( 10) .HRS( 10) .RR( 10) .PO( 10 
1).PT(10).PKSS(10).TDF(10) 

COMMON  /</  NML( 10) .POL( 10) .PTL( 10) .TOL( 10) .SRL( 10) .SVL( 10) .OZL( 10) 
I .RRL( 10) .HRSL( 10) .PKSL( 10) . ITITLE( 1 1 ) .TOLF( 10) .NLCHRL( 10) .TERlNO( I 
SO) 

COMMON  /PS/  NPENT.NPT( 10) .SALSZ( 10) .NLCHR( 10) .PLCHROY( 10) .ROS( 10) 
COMMON  /P3/  PENC.TR.SP.CORWD THIS) .DEL 

COMMON  /PL/  PSTART. ICHAFF.CHAFFS.CHAFFE. lECM.ECMS.ECME. IDECOt .uSTR 
I.OENO.CHEFF(IO) .ECMEFF(IO) 

COMMON  /PLH/  IPINT.0F( 10) .0ECPPEN( 10) .PENCDEC 
COMMON  /NHLN/  NSAM (66 . 1 0 . S) 

DIMENSION  IC0M(70) 

DATA  ( ( I COM (I ) . I ■ I . 5L ) -LHNPH I . 5HNPENT . LHPENC . SHTR . SHS? . SH  . SHHM . S 
I HVP . 3HROS . 3HNML . 7HEN0CASE . 6HEN0JCB . 6H  . 5HSALSZ . 7HPLCHROY . LH 

S .3H0EL.HHIL0W.5HALPHA.SH  .5H  .SH  .5HN01ST.6H  .SH  .hh 

3 1 SET . 7H  . 5H 1 SETL . SHASPEC  T . 3HPKH . 3HPKL . 5HT 1 TLE . 6HPSTART . 6H I CH 

HAFF . 6HCHAFFS . 6HCHAFFE . SHCHEFF . HH  I ECM . HHECMS . HHEC.ME . 6HECMEFF . 6H I DEC 
SOY . HHOSTR . HHOENO . 6HTER I NO . HHSRLO . 6H  . 7H0ECPPEN . 7HPENCOEC .SHIP 

6 1 NT . SHCF . THCORWOTH . 5HNSAML . SHNSAMH ) 

INITIALIZE  illegitimate  IDENTIFIER  FLAG 
IX«0 

EJECT  PAGE 
PRINT  33 

READ  IDENTIFIER  CARO 
I READ  3H.  NAME.X 
PRINT  35.  NAME.X 

BEGIN  CHECKING  FOR  IDENTIFIER  MATCH 
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IF  (NAME.EQ. ICOM( I ) ) NPHI-X 
IF  (NAME.NE.  lC0r*O)  1 GO  TO  2 
NPENT-X 

READ  36.  (NPTH).I-l.NPENT) 

PRINT  37.  (NPT( I ) . I«l .NPENT; 

2 continue 

IF  (NAHE.EQ. IC0M(3) ) PENC-X 
IF  (NAME.EQ. ICOM(H) » TR-X 
IF  (NAME.EQ. IC0M(5) ) SP-X 
IF  (NAME.EQ. IC0M(7) ) HM«X 
IF  (NAME.EQ. IC0M(8) ) VP-X 
IF  (NAME.NE. IC0M(9) > CO  TO  3 
READ  39.  (ROS(I).I-l.lO) 

PRINT  38.  (ROS(I) .1-1.10) 

3 CONTINUE 

IF  (NAME.NE.  ICOM(  10))  GO  TO  «+ 

REAO  36.  (NML(I).I-l.lO) 

PRINT  37.  (NML(I) .1-1.10) 

^ CONTINUE 

IF  (NAME.EQ. ICOM( 1 1 ) ) GO  TO  31 
IF  (NAME.EQ. ICOM( 12) ) CALL  EXIT 
IF  (NAME.NE.  1C0M(  m) ) GO  TO  5 
REAO  39.  (SALSZ(I).I-l.lO) 

PRINT  38.  (SALS2(I).I-1.10) 

5 continue 

IF  (NAME.NE. ICOM( 15) ) CO  TO  6 
REAO  39.  (PLCHROY(I).l-l.lO) 

PRINT  38.  (PLCHROYI I ) . I-l . 10) 

6 continue 

IF  (NAME.EQ. ICOM( 17))  CEL-X 
IF  (NAME.EQ. ICOM( 18) ) ILOU-X 
IF  (NAME.NE. ICOM( 19) ) GO  TO  7 
REAO  39.  (ALPHA! I >.1-1.10) 

PRINT  38.  (ALPHA! I ) . I-l . 10) 

7 continue 

IF  (NAME.EQ. IC0M(23) ) NOIST-X 
IF  (NAME.NE. IC0M(26) ) GO  TO  9 
C REAO  HIGH  altitude  SCENARIO  SAM  DATA 

REAO  HO.  (ISET(I).I-l.lO) 

PRINT  HI.  (IS£T(I).I-I.I0) 

PRINT  H2 
00  8 J-1 .10 

IF  (ISET(J).NE.I)  go  TO  8 
l-J 

REAO  H3.  I.NLCHR(I).PO(I).PT(I).TOF(1).TD(I).SR(1).SV(I).02(1).RR( 
11) .HRS(t) 

PRINT  H3.  I,NLCHR(I).PO(I),PT(I).TOF(I).TO(I),SR(1).SV(I).OZ(I) .RR 
Id)  .HRS(I) 

8 continue 
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9 CONTINUE 

IF  (NAME.NE. IC0M(28) ) GO  TO  11 
C REAO  LOW  altitude  SCENARIO  SAM  DATA 

READ  HO.  (ISETLI I) .1-1,10) 

PRINT  HI.  (ISETL(I) .1-1,10) 

PRINT  H2 
DO  10  J-1.10 

IF  (ISETL(J).NE.l)  GO  TO  10 
l-J 

READ  H3.  I .NLCHRLl I ) .POL( 1 ) .PTLf I ) .TDLF( I ) .TOLf • .SRL( I ) .SVL( 1 ) .DZ 
1L(1).RRLC1),HRSL(I) 

PRINT  H3.  I .NLCHRLl I ) .PDL( I ) .PTL( 1 ) .TOLF( I ) .TOL  ) .SRL( I ) .SVL( I ) .0 
1ZL(I).RRL(I).HRSL(I) 

10  CONTINUE 

11  continue 

IF  (NAME.NE. IC0M(29) ) GO  TO  13 
READ  39.  (ASPECT( I ) . I-l . 10) 

PRINT  38.  (ASPECT! I) .1-1.10) 

DO  12  I-l. 10 

ASPECT ( I ) -ASPECT ( 1 ) *6 . 283 1 853/ 360 . 0 

12  CONTINUE 

13  CONTINUE 

IF  (NAME.NE. ICOM(30) ) GO  TO  IH 
READ  39.  (PKSS! I). I-l. 10) 

PRINT  HH.  (PKSS(I).I-l.lO) 

iH  continue 

IF  (NAME.NE. 1C0M(31 ) ) GO  TO  15 
READ  39.  (PKSL(I).I-l.tO) 

PRINT  38,  (PKSL(I).l-l.lO) 

15  CONTINUE 

IF  (NAME.NE. IC0M(32) ) GO  TO  IS 
READ  H5.  (ITITLE(I).I-1.8) 

PRINT  H6.  (ITITLE(I).I-1.8) 

16  CONTINUE 

IF  (NAME.EQ, IC0M(33) ) PSTART-X 
IF  (NAME.EQ. IC0M(3H) ) ICHAFF-X 
IF  (NAME.EQ. IC0M(35) ) CHAFFS-X 
IF  (NAME.EQ. I COM (36) ) CHAFFE-X 
IF  (NAME.NE. IC0M(37) ) GO  TO  17 
REAO  39.  (CHEFF(I).I-I.IO) 

PRINT  30.  (CHEFFID.I-l.IO) 

17  CONTINUE 

IF  (NAME.EQ. IC0M(38) ) lECM-X 
IF  (NAME.EQ. 1C0M( 39) ) ECMS-X 
IF  (NAME.EQ. tCOM(HO) ) ECME-X 
IF  (NAME.NE. ICOM(Hl))  00  TO  18 
REAO  39.  (ECMEFF(I).I-l.lO) 

PRINT  38.  (ECMEFFin.I-l.tO) 

10  CONTINUE 

IF  (NAME.EQ. IC0M(H2) ) lOECOY-X 
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IF  (NAME.EQ. ICCM(43) ) DSTR-X 
IF  (NAME.EQ. 1C0M(^^) ) DENO-X 
IF  (NAME.NE. ICCM(H5) ) GO  TO  19 
READ  39.  (TERlNO(I).I-l.lO) 

PRINT  38.  (TERINOtn.l-l.lO) 

19  CONTINUE 

IF  (NAME.NE.  ICOMC+en  GO  TO  20 
READ  39.  (SRLD( I ) . 1-1 .10) 

PRINT  38.  (SRLD( I ) . I-l . 10) 

20  CONTINUE 

IF  (NAME.NE.  ICOMmS)  ) GO  TO  21 
READ  39.  (DECPPEN(I).I-l.lO) 

PRINT  38.  (OECPPENl 1 ) . I-I . 10) 

21  CONTINUE 

IF  (NAME.EQ.  ICCMm9)  ) PENCDEC-X 
IF  (NAME.EQ. ICOM(50) ) IPINT-X 
IF  (NAME.NE. IC0M(51 ) ) GO  TO  2^ 

IF  (X.LE.O.)  go  to  23 
DO  22  I-l. 10 

22  CF(I)-X 
GO  TO  2H 

23  CONTINUE 

READ  39.  (OF(  I ) .1-1.10 
PRINT  38.  (OF(I).I-l.lO) 

CONTINUE 

IF  (NAME.NE. IC0M(52) ) GO  TO  25 
READ  39.  CCRKOTHC 1 ) ,CORWOTH(2) 

PR  I NT  38 . CORWCTH ( 1 ) . CORWOTH ( 2 ) 

25  CONTINUE 

IF  (NAME.NE. 1CCM(53) .A. NAME. NE. ICOMtSH) ) GO  TO  29 
C READ  DEFENSE  SCENARIO  (250 

K-1 

IF  (NAME.EQ. IC0M(5H) ) K-2 

READ  36.  NNL 

PRINT  36.  NNL 

PRINT  32 

NNN-X 

DO  26  I-l  .66 

26  NSAM( I .NNN.K)-0 
DO  27  1-1, NNL 

READ  38.  NZCN.NSAM(NZON^l .NNN.K) 

27  PRINT  36.  NZON.NSAMlNZON^l .NNN.K) 

DO  28  1-2.66 

28  NSAM(  I .NNN.K)-NSAM(  I ,NNN  .K ) -^NSAMl  I - 1 .NNN  .K ) 

29  CONTINUE 

C DETERMINE  IF  IDENTIFIER  IS  LEGITIMATE 

DO  30  I-l .5^ 

IF  ( 1CCM( I ) .EQ.NAME)  go  to  1 

30  continue 

C IDENTIFIER  IS  NOT  LEGITIMATE 
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BBT  AVAIWBIE  COPY 


ooo 


IX-l 

PRINT  H7.  NAME 
CO  TO  I 

C EXIT  IF  AT  LEAST  ONE  ILLEGITIMATE  IDENTIFIER 

31  IF  (IX.EQ.l)  CALL  EXIT 
RETURN 


32  FORMAT  {13H  ZONE  NUMBER) 

33  FORMAT  (IHl) 

3V  FORMAT  (A7.3X.E20.8) 

35  FORMAT  ( IOX.A7.3X.F10.3) 

36  FORMAT  (1015) 

37  FORMAT  (10X.10I5) 

38  FORMAT  ( 1 OX. I OF  10. 2) 

39  FORMAT  (10F8.0) 

VO  FORMAT  (1011) 

VI  FORMAT  ( ICX.  lOI  I ) 

va  FORMAT  (7H  T NL/S.6X .2HP0 .5X .2HPT .5X . 3HTDF , 5X .2HT0 .6X . 2HSR . SX . 2HSV 
I . 5X . 2H0Z . VX . 2HRR . 8X . 3HHRS ) 

V3  FORMAT  ( 12. I3.3X.8F8.2.F8.3) 

VV  FORMAT  ( lOX. IOFIO.2) 

V5  FORMAT  (8A10) 

V6  FORMAT  (1QX.9A10) 
j V7  FORMAT  (2X.A7.21H  IS  NOT  IN  DICTIONARY) 

i END 


BEST  4V4IMBIE  COPY 
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